Global emergency and disaster transmission

ABSTRACT

Disclosed is a system and method for disaster warning recovery including a power modulator. A location is determined for an area affected by an emergency event. An emergency message is generated that corresponds to the emergency event. The emergency message is transmitted by via a guided surface wave. The guided surface wave is launched by a guided surface waveguide probe. The power modulator is coupled to the guided surface waveguide probe.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation of co-pending U.S. Non-provisionalPatent Application entitled “Global Emergency and DisasterTransmission,” which was filed on Sep. 8, 2015 and assigned applicationSer. No. 14/847,310, and was published on Mar. 9, 2017 as PublicationNumber US2017/0069189. This application is related to co-pending U.S.Non-provisional Patent Application entitled “Excitation and Use ofGuided Surface Wave Modes on Lossy Media,” which was filed on Mar. 7,2013 and assigned application Ser. No. 13/789,538, and was published onSep. 11, 2014 as Publication Number US2014/0252886 A1, and which isincorporated herein by reference in its entirety. This application isalso related to co-pending U.S. Non-provisional Patent Applicationentitled “Excitation and Use of Guided Surface Wave Modes on LossyMedia,” which was filed on Mar. 7, 2013 and assigned application Ser.No. 13/789,525, and was published on Sep. 11, 2014 as Publication NumberUS2014/0252865 A1, and which is incorporated herein by reference in itsentirety. This application is further related to co-pending U.S.Non-provisional patent application entitled “Excitation and Use ofGuided Surface Wave Modes on Lossy Media,” which was filed on Sep. 10,2014 and assigned application Ser. No. 14/483,089, and which isincorporated herein by reference in its entirety. This application isfurther related to co-pending U.S. Non-provisional patent applicationentitled “Excitation and Use of Guided Surface Waves,” which was filedon Jun. 2, 2015 and assigned application Ser. No. 14/728,507, and whichis incorporated herein by reference in its entirety. This application isfurther related to co-pending U.S. Non-provisional patent applicationentitled “Excitation and Use of Guided Surface Waves,” which was filedon Jun. 2, 2015 and assigned application Ser. No. 14/728,492, and whichis incorporated herein by reference in its entirety.

BACKGROUND

When disaster events occur, time is of the essence. Disaster events canoccur anywhere, including places without power. Numerous ways to detectthe occurrence of disaster events exist, such as monitoringseismographic data, weather patterns, and reports by witnesses to theevent, among other known methods. However, it can be difficult to notifypeople in the proximate area to the event, especially when the proximatearea is not near a power source. This problem may be exacerbated inthird-world countries or areas where power is not reliable. Further,even where power is readily and reliably available, when a disasterevent takes place, power in the affected area may not be reliable.

BRIEF DESCRIPTION OF THE DRAWINGS

Many aspects of the present disclosure can be better understood withreference to the following drawings. The components in the drawings arenot necessarily to scale, emphasis instead being placed upon clearlyillustrating the principles of the disclosure. Moreover, in thedrawings, like reference numerals designate corresponding partsthroughout the several views.

FIG. 1 is a chart that depicts field strength as a function of distancefor a guided electromagnetic field and a radiated electromagnetic field.

FIG. 2 is a drawing that illustrates a propagation interface with tworegions employed for transmission of a guided surface wave according tovarious embodiments of the present disclosure.

FIG. 3 is a drawing that illustrates a guided surface waveguide probedisposed with respect to a propagation interface of FIG. 2 according tovarious embodiments of the present disclosure.

FIG. 4 is a plot of an example of the magnitudes of close-in and far-outasymptotes of first order Hankel functions according to variousembodiments of the present disclosure.

FIGS. 5A and 5B are drawings that illustrate a complex angle ofincidence of an electric field synthesized by a guided surface waveguideprobe according to various embodiments of the present disclosure.

FIG. 6 is a graphical representation illustrating the effect ofelevation of a charge terminal on the location where the electric fieldof FIG. 5A intersects with the lossy conducting medium at a Brewsterangle according to various embodiments of the present disclosure.

FIG. 7 is a graphical representation of an example guided surfacewaveguide probe according to various embodiments of the presentdisclosure.

FIGS. 8A through 8C are graphical representations illustrating examplesof equivalent image plane models of the guided surface waveguide probeof FIGS. 3 and 7 according to various embodiments of the presentdisclosure.

FIGS. 9A and 9B are graphical representations illustrating examplemodels of a single-wire transmission line and a classic transmissionline of the equivalent image plane models of FIGS. 8B and 8C accordingto various embodiments of the present disclosure.

FIG. 10 is a flow chart illustrating an example of adjusting a guidedsurface waveguide probe of FIGS. 3 and 7 to launch a guided surface wavealong the surface of a lossy conducting medium according to variousembodiments of the present disclosure.

FIG. 11 is a plot illustrating an example of the relationship between awave tilt angle and the phase delay of a guided surface waveguide probeof FIGS. 3 and 7 according to various embodiments of the presentdisclosure.

FIG. 12 is a drawing that illustrates an example of a guided surfacewaveguide probe according to various embodiments of the presentdisclosure.

FIG. 13 is a graphical representation illustrating the incidence of asynthesized electric field at a complex Brewster angle to match theguided surface waveguide mode at the Hankel crossover distance accordingto various embodiments of the present disclosure.

FIG. 14 is a graphical representation of an example of a guided surfacewaveguide probe of FIG. 12 according to various embodiments of thepresent disclosure.

FIG. 15A includes plots of an example of the imaginary and real parts ofa phase delay (Φ_(U)) of a charge terminal T₁ of a guided surfacewaveguide probe according to various embodiments of the presentdisclosure.

FIG. 15B is a schematic diagram of the guided surface waveguide probe ofFIG. 14 according to various embodiments of the present disclosure.

FIG. 16 is a drawing that illustrates an example of a guided surfacewaveguide probe according to various embodiments of the presentdisclosure.

FIG. 17 is a graphical representation of an example of a guided surfacewaveguide probe of FIG. 16 according to various embodiments of thepresent disclosure.

FIGS. 18A through 18C depict examples of receiving structures that canbe employed to receive energy transmitted in the form of a guidedsurface wave launched by a guided surface waveguide probe according tothe various embodiments of the present disclosure.

FIG. 18D is a flow chart illustrating an example of adjusting areceiving structure according to various embodiments of the presentdisclosure.

FIG. 19 depicts an example of an additional receiving structure that canbe employed to receive energy transmitted in the form of a guidedsurface wave launched by a guided surface waveguide probe according tothe various embodiments of the present disclosure.

FIGS. 20A through 20E depict examples of various schematic symbolsaccording to various embodiments of the present disclosure.

FIG. 21 is an illustration of a networked environment according tovarious embodiments of the present disclosure.

FIG. 22 depicts a disaster warning device according to variousembodiments of the present disclosure.

FIG. 23 depicts a disaster warning device according to variousembodiments of the present disclosure.

FIG. 24 is a flowchart illustrating one example of functionalityimplemented as portions of an emergency application executed in acomputing environment of FIG. 21 according to various embodiments of thepresent disclosure.

FIG. 25 is a flowchart illustrating one example of functionalityimplemented as portions of an emergency application executed in acomputing environment of FIG. 21 according to various embodiments of thepresent disclosure.

FIG. 26 is a flowchart illustrating one example of functionalityimplemented as portions of a disaster application 537 of FIGS. 22 and 23according to various embodiments of the present disclosure.

FIG. 27 is a flowchart illustrating one example of functionalityimplemented as portions of a disaster application 537 of FIGS. 22 and 23according to various embodiments of the present disclosure.

FIG. 28 is a flowchart illustrating one example of functionalityimplemented as portions of an emergency application executed in acomputing environment of FIG. 21 according to various embodiments of thepresent disclosure.

FIG. 29 is a schematic block diagram that provides one exampleillustration of a computing environment employed in the networkedenvironment of FIG. 21 according to various embodiments of the presentdisclosure.

DETAILED DESCRIPTION

To begin, some terminology shall be established to provide clarity inthe discussion of concepts to follow. First, as contemplated herein, aformal distinction is drawn between radiated electromagnetic fields andguided electromagnetic fields.

As contemplated herein, a radiated electromagnetic field compriseselectromagnetic energy that is emitted from a source structure in theform of waves that are not bound to a waveguide. For example, a radiatedelectromagnetic field is generally a field that leaves an electricstructure such as an antenna and propagates through the atmosphere orother medium and is not bound to any waveguide structure. Once radiatedelectromagnetic waves leave an electric structure such as an antenna,they continue to propagate in the medium of propagation (such as air)independent of their source until they dissipate regardless of whetherthe source continues to operate. Once electromagnetic waves areradiated, they are not recoverable unless intercepted, and, if notintercepted, the energy inherent in the radiated electromagnetic wavesis lost forever. Electrical structures such as antennas are designed toradiate electromagnetic fields by maximizing the ratio of the radiationresistance to the structure loss resistance. Radiated energy spreads outin space and is lost regardless of whether a receiver is present. Theenergy density of the radiated fields is a function of distance due togeometric spreading. Accordingly, the term “radiate” in all its forms asused herein refers to this form of electromagnetic propagation.

A guided electromagnetic field is a propagating electromagnetic wavewhose energy is concentrated within or near boundaries between mediahaving different electromagnetic properties. In this sense, a guidedelectromagnetic field is one that is bound to a waveguide and may becharacterized as being conveyed by the current flowing in the waveguide.If there is no load to receive and/or dissipate the energy conveyed in aguided electromagnetic wave, then no energy is lost except for thatdissipated in the conductivity of the guiding medium. Stated anotherway, if there is no load for a guided electromagnetic wave, then noenergy is consumed. Thus, a generator or other source generating aguided electromagnetic field does not deliver real power unless aresistive load is present. To this end, such a generator or other sourceessentially runs idle until a load is presented. This is akin to runninga generator to generate a 60 Hertz electromagnetic wave that istransmitted over power lines where there is no electrical load. Itshould be noted that a guided electromagnetic field or wave is theequivalent to what is termed a “transmission line mode.” This contrastswith radiated electromagnetic waves in which real power is supplied atall times in order to generate radiated waves. Unlike radiatedelectromagnetic waves, guided electromagnetic energy does not continueto propagate along a finite length waveguide after the energy source isturned off. Accordingly, the term “guide” in all its forms as usedherein refers to this transmission mode of electromagnetic propagation.

Referring now to FIG. 1, shown is a graph 100 of field strength indecibels (dB) above an arbitrary reference in volts per meter as afunction of distance in kilometers on a log-dB plot to furtherillustrate the distinction between radiated and guided electromagneticfields. The graph 100 of FIG. 1 depicts a guided field strength curve103 that shows the field strength of a guided electromagnetic field as afunction of distance. This guided field strength curve 103 isessentially the same as a transmission line mode. Also, the graph 100 ofFIG. 1 depicts a radiated field strength curve 106 that shows the fieldstrength of a radiated electromagnetic field as a function of distance.

Of interest are the shapes of the curves 103 and 106 for guided wave andfor radiation propagation, respectively. The radiated field strengthcurve 106 falls off geometrically (1/d, where d is distance), which isdepicted as a straight line on the log-log scale. The guided fieldstrength curve 103, on the other hand, has a characteristic exponentialdecay of e^(−ad)/√{square root over (d)} and exhibits a distinctive knee109 on the log-log scale.

The guided field strength curve 103 and the radiated field strengthcurve 106 intersect at point 112, which occurs at a crossing distance.At distances less than the crossing distance at intersection point 112,the field strength of a guided electromagnetic field is significantlygreater at most locations than the field strength of a radiatedelectromagnetic field. At distances greater than the crossing distance,the opposite is true. Thus, the guided and radiated field strengthcurves 103 and 106 further illustrate the fundamental propagationdifference between guided and radiated electromagnetic fields. For aninformal discussion of the difference between guided and radiatedelectromagnetic fields, reference is made to Milligan, T., ModernAntenna Design, McGraw-Hill, 1^(st) Edition, 1985, pp. 8-9, which isincorporated herein by reference in its entirety.

The distinction between radiated and guided electromagnetic waves, madeabove, is readily expressed formally and placed on a rigorous basis.That two such diverse solutions could emerge from one and the samelinear partial differential equation, the wave equation, analyticallyfollows from the boundary conditions imposed on the problem. The Greenfunction for the wave equation, itself, contains the distinction betweenthe nature of radiation and guided waves.

In empty space, the wave equation is a differential operator whoseeigenfunctions possess a continuous spectrum of eigenvalues on thecomplex wave-number plane. This transverse electro-magnetic (TEM) fieldis called the radiation field, and those propagating fields are called“Hertzian waves.” However, in the presence of a conducting boundary, thewave equation plus boundary conditions mathematically lead to a spectralrepresentation of wave-numbers composed of a continuous spectrum plus asum of discrete spectra. To this end, reference is made to Sommerfeld,A., “Uber die Ausbreitung der Wellen in der Drahtlosen Telegraphie,”Annalen der Physik, Vol. 28, 1909, pp. 665-736. Also see Sommerfeld, A.,“Problems of Radio,” published as Chapter 6 in Partial DifferentialEquations in Physics-Lectures on Theoretical Physics: Volume VI,Academic Press, 1949, pp. 236-289, 295-296; Collin, R. E., “HertzianDipole Radiating Over a Lossy Earth or Sea: Some Early and Late 20^(th)Century Controversies,” IEEE Antennas and Propagation Maqazine, Vol. 46,No. 2, April 2004, pp. 64-79; and Reich, H. J., Ordnung, P. F, Krauss,H. L., and Skalnik, J. G., Microwave Theory and Techniques, VanNostrand, 1953, pp. 291-293, each of these references being incorporatedherein by reference in its entirety.

The terms “ground wave” and “surface wave” identify two distinctlydifferent physical propagation phenomena. A surface wave arisesanalytically from a distinct pole yielding a discrete component in theplane wave spectrum. See, e.g., “The Excitation of Plane Surface Waves”by Cullen, A. L., (Proceedinqs of the IEE (British), Vol. 101, Part IV,August 1954, pp. 225-235). In this context, a surface wave is consideredto be a guided surface wave. The surface wave (in the Zenneck-Sommerfeldguided wave sense) is, physically and mathematically, not the same asthe ground wave (in the Weyl-Norton-FCC sense) that is now so familiarfrom radio broadcasting. These two propagation mechanisms arise from theexcitation of different types of eigenvalue spectra (continuum ordiscrete) on the complex plane. The field strength of the guided surfacewave decays exponentially with distance as illustrated by curve 103 ofFIG. 1 (much like propagation in a lossy waveguide) and resemblespropagation in a radial transmission line, as opposed to the classicalHertzian radiation of the ground wave, which propagates spherically,possesses a continuum of eigenvalues, falls off geometrically asillustrated by curve 106 of FIG. 1, and results from branch-cutintegrals. As experimentally demonstrated by C. R. Burrows in “TheSurface Wave in Radio Propagation over Plane Earth” (Proceedinqs of theIRE, Vol. 25, No. 2, February, 1937, pp. 219-229) and “The Surface Wavein Radio Transmission” (Bell Laboratories Record, Vol. 15, June 1937,pp. 321-324), vertical antennas radiate ground waves but do not launchguided surface waves.

To summarize the above, first, the continuous part of the wave-numbereigenvalue spectrum, corresponding to branch-cut integrals, produces theradiation field, and second, the discrete spectra, and correspondingresidue sum arising from the poles enclosed by the contour ofintegration, result in non-TEM traveling surface waves that areexponentially damped in the direction transverse to the propagation.Such surface waves are guided transmission line modes. For furtherexplanation, reference is made to Friedman, B., Principles andTechniques of Applied Mathematics, Wiley, 1956, pp. pp. 214, 283-286,290, 298-300.

In free space, antennas excite the continuum eigenvalues of the waveequation, which is a radiation field, where the outwardly propagating RFenergy with E_(z) and H_(ϕ) in-phase is lost forever. On the other hand,waveguide probes excite discrete eigenvalues, which results intransmission line propagation. See Collin, R. E., Field Theory of GuidedWaves, McGraw-Hill, 1960, pp. 453, 474-477. While such theoreticalanalyses have held out the hypothetical possibility of launching opensurface guided waves over planar or spherical surfaces of lossy,homogeneous media, for more than a century no known structures in theengineering arts have existed for accomplishing this with any practicalefficiency. Unfortunately, since it emerged in the early 1900's, thetheoretical analysis set forth above has essentially remained a theoryand there have been no known structures for practically accomplishingthe launching of open surface guided waves over planar or sphericalsurfaces of lossy, homogeneous media.

According to the various embodiments of the present disclosure, variousguided surface waveguide probes are described that are configured toexcite electric fields that couple into a guided surface waveguide modealong the surface of a lossy conducting medium. Such guidedelectromagnetic fields are substantially mode-matched in magnitude andphase to a guided surface wave mode on the surface of the lossyconducting medium. Such a guided surface wave mode can also be termed aZenneck waveguide mode. By virtue of the fact that the resultant fieldsexcited by the guided surface waveguide probes described herein aresubstantially mode-matched to a guided surface waveguide mode on thesurface of the lossy conducting medium, a guided electromagnetic fieldin the form of a guided surface wave is launched along the surface ofthe lossy conducting medium. According to one embodiment, the lossyconducting medium comprises a terrestrial medium such as the Earth.

Referring to FIG. 2, shown is a propagation interface that provides foran examination of the boundary value solutions to Maxwell's equationsderived in 1907 by Jonathan Zenneck as set forth in his paper Zenneck,J., “On the Propagation of Plane Electromagnetic Waves Along a FlatConducting Surface and their Relation to Wireless Telegraphy,” Annalender Physik, Serial 4, Vol. 23, Sep. 20, 1907, pp. 846-866. FIG. 2depicts cylindrical coordinates for radially propagating waves along theinterface between a lossy conducting medium specified as Region 1 and aninsulator specified as Region 2. Region 1 can comprise, for example, anylossy conducting medium. In one example, such a lossy conducting mediumcan comprise a terrestrial medium such as the Earth or other medium.Region 2 is a second medium that shares a boundary interface with Region1 and has different constitutive parameters relative to Region 1. Region2 can comprise, for example, any insulator such as the atmosphere orother medium. The reflection coefficient for such a boundary interfacegoes to zero only for incidence at a complex Brewster angle. SeeStratton, J. A., Electromaqnetic Theory, McGraw-Hill, 1941, p. 516.

According to various embodiments, the present disclosure sets forthvarious guided surface waveguide probes that generate electromagneticfields that are substantially mode-matched to a guided surface waveguidemode on the surface of the lossy conducting medium comprising Region 1.According to various embodiments, such electromagnetic fieldssubstantially synthesize a wave front incident at a complex Brewsterangle of the lossy conducting medium that can result in zero reflection.

To explain further, in Region 2, where an e^(jωt) field variation isassumed and where ρ≠0 and z≥0 (with z being the vertical coordinatenormal to the surface of Region 1, and ρ being the radial dimension incylindrical coordinates), Zenneck's closed-form exact solution ofMaxwell's equations satisfying the boundary conditions along theinterface are expressed by the following electric field and magneticfield components:

$\begin{matrix}{{H_{2\varphi} = {A\; e^{{- u_{2}}z}{H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}}},} & (1) \\{{E_{2\rho} = {{A\left( \frac{u_{2}}{j\; \omega \; ɛ_{o}} \right)}e^{{- u_{2}}z}{H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}}},{and}} & (2) \\{E_{2z} = {{A\left( \frac{- \gamma}{\omega \; ɛ_{o}} \right)}e^{{- u_{2}}z}{{H_{0}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}.}}} & (3)\end{matrix}$

In Region 1, where the e^(jωt) field variation is assumed and where ρ≠0and z≤0, Zenneck's closed-form exact solution of Maxwell's equationssatisfying the boundary conditions along the interface is expressed bythe following electric field and magnetic field components:

$\begin{matrix}{{H_{1\varphi} = {A\; e^{u_{1}z}{H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}}},} & (4) \\{{E_{1\rho} = {{A\left( \frac{- u_{1}}{\sigma_{1} + {j\; \omega \; ɛ_{1}}} \right)}e^{u_{1}z}{H_{1}^{(2)}\left( {{- j}\; {\gamma\rho}} \right)}}},{and}} & (5) \\{E_{1z} = {{A\left( \frac{{- j}\; \gamma}{\sigma_{1} + {j\; \omega \; ɛ_{1}}} \right)}e^{u_{1}z}{{H_{1}^{(2)}\left( {{- j}\; {\gamma\rho}} \right)}.}}} & (6)\end{matrix}$

In these expressions, z is the vertical coordinate normal to the surfaceof Region 1 and p is the radial coordinate, H_(n) ⁽²⁾(−jγρ) is a complexargument Hankel function of the second kind and order n, u₁ is thepropagation constant in the positive vertical (z) direction in Region 1,u₂ is the propagation constant in the vertical (z) direction in Region2, σ₁ is the conductivity of Region 1, ω is equal to 2πf, where f is afrequency of excitation, ε₀ is the permittivity of free space, ε₁ is thepermittivity of Region 1, A is a source constant imposed by the source,and γ is a surface wave radial propagation constant.

The propagation constants in the ±z directions are determined byseparating the wave equation above and below the interface betweenRegions 1 and 2, and imposing the boundary conditions. This exercisegives, in Region 2,

$\begin{matrix}{u_{2} = \frac{- {jk}_{o}}{\sqrt{1 + \left( {ɛ_{r} - {jx}} \right)}}} & (7)\end{matrix}$

and gives, in Region 1,

u ₁ =−u ₂(ε_(r) −jx).  (8)

The radial propagation constant γ is given by

$\begin{matrix}{{\gamma = {{j\sqrt{k_{o}^{2} + u_{2}^{2}}} = {j\frac{\; {k_{o}n}}{\sqrt{1 + n^{2}}}}}},} & (9)\end{matrix}$

which is a complex expression where n is the complex index of refractiongiven by

n=√{square root over (ε_(r) −jx)}.  (10)

In all of the above Equations,

$\begin{matrix}{{x = \frac{\sigma_{1}}{\omega \; ɛ_{o}}},{and}} & (11) \\{{k_{o} = {{\omega \sqrt{\mu_{o}ɛ_{o}}} = \frac{\lambda_{o}}{2\pi}}},} & (12)\end{matrix}$

where ε_(r) comprises the relative permittivity of Region 1, σ₁ is theconductivity of Region 1, ε_(o) is the permittivity of free space, andμ_(o) comprises the permeability of free space. Thus, the generatedsurface wave propagates parallel to the interface and exponentiallydecays vertical to it. This is known as evanescence.

Thus, Equations (1)-(3) can be considered to be acylindrically-symmetric, radially-propagating waveguide mode. SeeBarlow, H. M., and Brown, J., Radio Surface Waves, Oxford UniversityPress, 1962, pp. 10-12, 29-33. The present disclosure details structuresthat excite this “open boundary” waveguide mode. Specifically, accordingto various embodiments, a guided surface waveguide probe is providedwith a charge terminal of appropriate size that is fed with voltageand/or current and is positioned relative to the boundary interfacebetween Region 2 and Region 1. This may be better understood withreference to FIG. 3, which shows an example of a guided surfacewaveguide probe 200 a that includes a charge terminal T₁ elevated abovea lossy conducting medium 203 (e.g., the Earth) along a vertical axis zthat is normal to a plane presented by the lossy conducting medium 203.The lossy conducting medium 203 makes up Region 1, and a second medium206 makes up Region 2 and shares a boundary interface with the lossyconducting medium 203.

According to one embodiment, the lossy conducting medium 203 cancomprise a terrestrial medium such as the planet Earth. To this end,such a terrestrial medium comprises all structures or formationsincluded thereon whether natural or man-made. For example, such aterrestrial medium can comprise natural elements such as rock, soil,sand, fresh water, sea water, trees, vegetation, and all other naturalelements that make up our planet. In addition, such a terrestrial mediumcan comprise man-made elements such as concrete, asphalt, buildingmaterials, and other man-made materials. In other embodiments, the lossyconducting medium 203 can comprise some medium other than the Earth,whether naturally occurring or man-made. In other embodiments, the lossyconducting medium 203 can comprise other media such as man-made surfacesand structures such as automobiles, aircraft, man-made materials (suchas plywood, plastic sheeting, or other materials) or other media.

In the case where the lossy conducting medium 203 comprises aterrestrial medium or Earth, the second medium 206 can comprise theatmosphere above the ground. As such, the atmosphere can be termed an“atmospheric medium” that comprises air and other elements that make upthe atmosphere of the Earth. In addition, it is possible that the secondmedium 206 can comprise other media relative to the lossy conductingmedium 203.

The guided surface waveguide probe 200 a includes a feed network 209that couples an excitation source 212 to the charge terminal T₁ via,e.g., a vertical feed line conductor. According to various embodiments,a charge Q₁ is imposed on the charge terminal T₁ to synthesize anelectric field based upon the voltage applied to terminal T₁ at anygiven instant. Depending on the angle of incidence (θ_(i)) of theelectric field (E), it is possible to substantially mode-match theelectric field to a guided surface waveguide mode on the surface of thelossy conducting medium 203 comprising Region 1.

By considering the Zenneck closed-form solutions of Equations (1)-(6),the Leontovich impedance boundary condition between Region 1 and Region2 can be stated as

{circumflex over (z)}×{right arrow over (H)} ₂(ρ,φ,0)={right arrow over(J)} _(s),  (13)

where {circumflex over (z)} is a unit normal in the positive vertical(+z) direction and {right arrow over (H)}₂ is the magnetic fieldstrength in Region 2 expressed by Equation (1) above. Equation (13)implies that the electric and magnetic fields specified in Equations(1)-(3) may result in a radial surface current density along theboundary interface, where the radial surface current density can bespecified by

J _(ρ)(ρ′)=−AH ₁ ⁽²⁾(−jγρ′)  (14)

where A is a constant. Further, it should be noted that close-in to theguided surface waveguide probe 200 (for ρ<<λ), Equation (14) above hasthe behavior

$\begin{matrix}{{J_{close}\left( \rho^{\prime} \right)} = {\frac{- {A\left( {j\; 2} \right)}}{\pi \left( {{- j}\; \gamma \; \rho^{\prime}} \right)} = {{- H_{\varphi}} = {- {\frac{I_{o}}{2\; \pi \; \rho^{\prime}}.}}}}} & (15)\end{matrix}$

The negative sign means that when source current (I_(o)) flowsvertically upward as illustrated in FIG. 3, the “close-in” groundcurrent flows radially inward. By field matching on H_(ϕ) “close-in,” itcan be determined that

$\begin{matrix}{A = {{- \frac{I_{o}\gamma}{4}} = {- \frac{\omega \; q_{1}\gamma}{4}}}} & (16)\end{matrix}$

where q₁=C₁V₁, in Equations (1)-(6) and (14). Therefore, the radialsurface current density of Equation (14) can be restated as

$\begin{matrix}{{J_{\rho}\left( \rho^{\prime} \right)} = {\frac{I_{o}\gamma}{4}{{H_{1}^{(2)}\left( {{- j}\; {\gamma\rho}^{\prime}} \right)}.}}} & (17)\end{matrix}$

The fields expressed by Equations (1)-(6) and (17) have the nature of atransmission line mode bound to a lossy interface, not radiation fieldsthat are associated with groundwave propagation. See Barlow, H. M. andBrown, J., Radio Surface Waves, Oxford University Press, 1962, pp. 1-5.

At this point, a review of the nature of the Hankel functions used inEquations (1)-(6) and (17) is provided for these solutions of the waveequation. One might observe that the Hankel functions of the first andsecond kind and order n are defined as complex combinations of thestandard Bessel functions of the first and second kinds

H _(n) ⁽¹⁾(x)=J _(n)(x)+jN _(n)(x), and  (18)

H _(n) ⁽²⁾(x)=J _(n)(x)−jN _(n)(x),  (19)

These functions represent cylindrical waves propagating radially inward(H_(n) ⁽¹⁾) and outward (H_(n) ⁽²⁾), respectively. The definition isanalogous to the relationship e^(±jx)=cos x±j sin x. See, for example,Harrington, R. F., Time-Harmonic Fields, McGraw-Hill, 1961, pp. 460-463.

That H_(n) ⁽²⁾(k_(ρ)ρ) is an outgoing wave can be recognized from itslarge argument asymptotic behavior that is obtained directly from theseries definitions of J_(n)(x) and N_(n)(x). Far-out from the guidedsurface waveguide probe:

$\begin{matrix}{{{{H_{n}^{(2)}(x)}\underset{x\rightarrow\infty}{\rightarrow}{\sqrt{\frac{2j}{\pi \; x}}j^{n}e^{{- j}\; x}}} = {\sqrt{\frac{2}{\pi \; x}}j^{n}e^{- {j{({x - \frac{\pi}{4}})}}}}},} & \left( {20a} \right)\end{matrix}$

which, when multiplied by e^(Jωt), is an outward propagating cylindricalwave of the form e^(j(ωt-kρ)) with a 1/√{square root over (ρ)} spatialvariation. The first order (n=1) solution can be determined fromEquation (20a) to be

$\begin{matrix}{{{H_{1}^{(2)}(x)}\underset{x\rightarrow\infty}{\rightarrow}{j\sqrt{\frac{2j}{\pi \; x}}e^{{- j}\; x}}} = {\sqrt{\frac{2}{\pi \; x}}{e^{- {j{({x - \frac{\pi}{2} - \frac{\pi}{4}})}}}.}}} & \left( {20b} \right)\end{matrix}$

Close-in to the guided surface waveguide probe (for ρ<<λ), the Hankelfunction of first order and the second kind behaves as

$\begin{matrix}{{H_{1}^{(2)}(x)}\underset{x\rightarrow\infty}{\rightarrow}{\frac{2j}{\pi \; x}.}} & (21)\end{matrix}$

Note that these asymptotic expressions are complex quantities. When x isa real quantity, Equations (20b) and (21) differ in phase by √{squareroot over (j)}, which corresponds to an extra phase advance or “phaseboost” of 45° or, equivalently, λ/8. The close-in and far-out asymptotesof the first order Hankel function of the second kind have a Hankel“crossover” or transition point where they are of equal magnitude at adistance of ρ=R_(x).

Thus, beyond the Hankel crossover point the “far out” representationpredominates over the “close-in” representation of the Hankel function.The distance to the Hankel crossover point (or Hankel crossoverdistance) can be found by equating Equations (20b) and (21) for −jγρ,and solving for R_(x). With x=σ/ωε_(o), it can be seen that the far-outand close-in Hankel function asymptotes are frequency dependent, withthe Hankel crossover point moving out as the frequency is lowered. Itshould also be noted that the Hankel function asymptotes may also varyas the conductivity (σ) of the lossy conducting medium changes. Forexample, the conductivity of the soil can vary with changes in weatherconditions.

Referring to FIG. 4, shown is an example of a plot of the magnitudes ofthe first order Hankel functions of Equations (20b) and (21) for aRegion 1 conductivity of σ=0.010 mhos/m and relative permittivityε_(r)=15, at an operating frequency of 1850 kHz. Curve 115 is themagnitude of the far-out asymptote of Equation (20b) and curve 118 isthe magnitude of the close-in asymptote of Equation (21), with theHankel crossover point 121 occurring at a distance of R_(x)=54 feet.While the magnitudes are equal, a phase offset exists between the twoasymptotes at the Hankel crossover point 121. It can also be seen thatthe Hankel crossover distance is much less than a wavelength of theoperation frequency.

Considering the electric field components given by Equations (2) and (3)of the Zenneck closed-form solution in Region 2, it can be seen that theratio of ε_(z) and ε_(p) asymptotically passes to

$\begin{matrix}{{\frac{E_{z}}{E_{\rho}} = {{{\left( \frac{{- j}\; \gamma}{u_{2}} \right)\frac{H_{0}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}{H_{1}^{(2)}\left( {{- j}\; {\gamma\rho}} \right)}}\underset{\rho\rightarrow\infty}{\rightarrow}\sqrt{ɛ_{r} - {j\; \frac{\sigma}{{\omega ɛ}_{o}}}}} = {n = {\tan \; \theta_{i}}}}},} & (22)\end{matrix}$

where n is the complex index of refraction of Equation (10) and θ_(i) isthe angle of incidence of the electric field. In addition, the verticalcomponent of the mode-matched electric field of Equation (3)asymptotically passes to

$\begin{matrix}{{E_{2z}\underset{\rho\rightarrow\infty}{\rightarrow}{\left( \frac{q_{free}}{ɛ_{o}} \right)\sqrt{\frac{\gamma^{3}}{8\pi}}e^{{- u_{2}}z}\frac{e^{- {j{({{\gamma\rho} - {\pi/4}})}}}}{\sqrt{\rho}}}},} & (23)\end{matrix}$

which is linearly proportional to free charge on the isolated componentof the elevated charge terminal's capacitance at the terminal voltage,q_(free)=C_(free)× V_(T).

For example, the height H₁ of the elevated charge terminal T₁ in FIG. 3affects the amount of free charge on the charge terminal T₁. When thecharge terminal T₁ is near the ground plane of Region 1, most of thecharge Q₁ on the terminal is “bound.” As the charge terminal T₁ iselevated, the bound charge is lessened until the charge terminal T₁reaches a height at which substantially all of the isolated charge isfree.

The advantage of an increased capacitive elevation for the chargeterminal T₁ is that the charge on the elevated charge terminal T₁ isfurther removed from the ground plane, resulting in an increased amountof free charge q_(free) to couple energy into the guided surfacewaveguide mode. As the charge terminal T₁ is moved away from the groundplane, the charge distribution becomes more uniformly distributed aboutthe surface of the terminal. The amount of free charge is related to theself-capacitance of the charge terminal T₁.

For example, the capacitance of a spherical terminal can be expressed asa function of physical height above the ground plane. The capacitance ofa sphere at a physical height of h above a perfect ground is given by

C _(elevated sphere)=4πε_(o) a(1+M+M ² +M ³+2M ⁴+3M ⁵+ . . . ),  (24)

where the diameter of the sphere is 2a, and where M=a/2h with h beingthe height of the spherical terminal. As can be seen, an increase in theterminal height h reduces the capacitance C of the charge terminal. Itcan be shown that for elevations of the charge terminal T₁ that are at aheight of about four times the diameter (4D=8a) or greater, the chargedistribution is approximately uniform about the spherical terminal,which can improve the coupling into the guided surface waveguide mode.

In the case of a sufficiently isolated terminal, the self-capacitance ofa conductive sphere can be approximated by C=4πε_(o)a, where a is theradius of the sphere in meters, and the self-capacitance of a disk canbe approximated by C=8ε_(o)a, where a is the radius of the disk inmeters. The charge terminal T₁ can include any shape such as a sphere, adisk, a cylinder, a cone, a torus, a hood, one or more rings, or anyother randomized shape or combination of shapes. An equivalent sphericaldiameter can be determined and used for positioning of the chargeterminal T₁.

This may be further understood with reference to the example of FIG. 3,where the charge terminal T₁ is elevated at a physical height ofh_(p)=H₁ above the lossy conducting medium 203. To reduce the effects ofthe “bound” charge, the charge terminal T₁ can be positioned at aphysical height that is at least four times the spherical diameter (orequivalent spherical diameter) of the charge terminal T₁ to reduce thebounded charge effects.

Referring next to FIG. 5A, shown is a ray optics interpretation of theelectric field produced by the elevated charge Q₁ on charge terminal T₁of FIG. 3. As in optics, minimizing the reflection of the incidentelectric field can improve and/or maximize the energy coupled into theguided surface waveguide mode of the lossy conducting medium 203. For anelectric field (E_(∥)) that is polarized parallel to the plane ofincidence (not the boundary interface), the amount of reflection of theincident electric field may be determined using the Fresnel reflectioncoefficient, which can be expressed as

$\begin{matrix}{{{\Gamma_{||}\left( \theta_{i} \right)} = {\frac{E_{||{,R}}}{E_{||{,i}}} = \frac{{\sqrt{\left. {ɛ_{r} - {jx}} \right) - {\sin^{2}\theta_{i}}} - {\left( {ɛ_{r} - {jx}} \right)\cos \; \theta_{i}}}\;}{\sqrt{\left( {ɛ_{r} - {jx}} \right) - {\sin^{2}\theta_{i}}} + {\left( {ɛ_{r} - {jx}} \right)\cos \; \theta_{i}}}}},} & (25)\end{matrix}$

where θ_(i) is the conventional angle of incidence measured with respectto the surface normal.

In the example of FIG. 5A, the ray optic interpretation shows theincident field polarized parallel to the plane of incidence having anangle of incidence of θ_(i) which is measured with respect to thesurface normal ({circumflex over (z)}). There will be no reflection ofthe incident electric field when Γ_(∥)(θ_(i))=0 and thus the incidentelectric field will be completely coupled into a guided surfacewaveguide mode along the surface of the lossy conducting medium 203. Itcan be seen that the numerator of Equation (25) goes to zero when theangle of incidence is

θ_(i)=arctan(√{square root over (ε_(r) −jx)})=θ_(i,B),  (26)

where x=σ/ωε_(o). This complex angle of incidence (θ_(i,B)) is referredto as the Brewster angle. Referring back to Equation (22), it can beseen that the same complex Brewster angle (θ_(i,B)) relationship ispresent in both Equations (22) and (26).

As illustrated in FIG. 5A, the electric field vector E can be depictedas an incoming non-uniform plane wave, polarized parallel to the planeof incidence. The electric field vector E can be created fromindependent horizontal and vertical components as

{right arrow over (E)}(θ_(i))=E _(ρ) {circumflex over (ρ)}+E _(z){circumflex over (z)}.  (27)

Geometrically, the illustration in FIG. 5A suggests that the electricfield vector E can be given by

E _(ρ)(ρ,z)=E(ρ,z)cos θ_(i), and  (28a)

E _(z)(ρ,z)=E(ρ,z)cos(π/2−θ_(i))=E(ρ,z)sin θ_(i),  (28b)

which means that the field ratio is

$\begin{matrix}{\frac{E_{\rho}}{E_{z}} = {\frac{1}{\tan \; \theta_{i}} = {\tan \; {\psi_{i}.}}}} & (29)\end{matrix}$

A generalized parameter W, called “wave tilt,” is noted herein as theratio of the horizontal electric field component to the verticalelectric field component given by

$\begin{matrix}{{W = {\frac{E_{\rho}}{E_{z}} = {{W}e^{j\; \Psi}}}},{or}} & \left( {30a} \right) \\{{\frac{1}{W} = {\frac{E_{z}}{E_{\rho}} = {{\tan \; \theta_{i}} = {\frac{1}{W}e^{{- j}\; \Psi}}}}},} & \left( {30b} \right)\end{matrix}$

which is complex and has both magnitude and phase. For anelectromagnetic wave in Region 2, the wave tilt angle (Ψ) is equal tothe angle between the normal of the wave-front at the boundary interfacewith Region 1 and the tangent to the boundary interface. This may beeasier to see in FIG. 5B, which illustrates equi-phase surfaces of anelectromagnetic wave and their normals for a radial cylindrical guidedsurface wave. At the boundary interface (z=0) with a perfect conductor,the wave-front normal is parallel to the tangent of the boundaryinterface, resulting in W=0. However, in the case of a lossy dielectric,a wave tilt W exists because the wave-front normal is not parallel withthe tangent of the boundary interface at z=0.

Applying Equation (30b) to a guided surface wave gives

$\begin{matrix}{{\tan \; \theta_{i,B}} = {\frac{E_{z}}{E_{\rho}} = {\frac{u_{2}}{\gamma} = {\sqrt{ɛ_{r} - {jx}} = {n = {\frac{1}{W} = {\frac{1}{W}{e^{{- j}\; \Psi}.}}}}}}}} & (31)\end{matrix}$

With the angle of incidence equal to the complex Brewster angle(θ_(i,B)), the Fresnel reflection coefficient of Equation (25) vanishes,as shown by

$\begin{matrix}{{\Gamma_{||}\left( \theta_{i,B} \right)} = {\left. \frac{\sqrt{\left( {ɛ_{r} - {jx}} \right) - {\sin^{2}\theta_{i}}} - {\left( {ɛ_{r} - {jx}} \right)\cos \; \theta_{i}}}{\sqrt{\left( {ɛ_{r} - {jx}} \right) - {\sin^{2}\theta_{i}}} + {\left( {ɛ_{r} - {jx}} \right)\cos \; \theta_{i}}} \right|_{\theta_{i} = \theta_{i,B}} = 0.}} & (32)\end{matrix}$

By adjusting the complex field ratio of Equation (22), an incident fieldcan be synthesized to be incident at a complex angle at which thereflection is reduced or eliminated. Establishing this ratio asn=√{square root over (ε_(r)−jx)} results in the synthesized electricfield being incident at the complex Brewster angle, making thereflections vanish.

The concept of an electrical effective height can provide furtherinsight into synthesizing an electric field with a complex angle ofincidence with a guided surface waveguide probe 200. The electricaleffective height (h_(eff)) has been defined as

$\begin{matrix}{h_{eff} = {\frac{1}{l_{0}}{\int_{0}^{h_{p}}{{I(z)}{dz}}}}} & (33)\end{matrix}$

for a monopole with a physical height (or length) of h_(p). Since theexpression depends upon the magnitude and phase of the sourcedistribution along the structure, the effective height (or length) iscomplex in general. The integration of the distributed current I(z) ofthe structure is performed over the physical height of the structure(h_(p)), and normalized to the ground current (I₀) flowing upwardthrough the base (or input) of the structure. The distributed currentalong the structure can be expressed by

I(z)=I _(C) cos(β₀ z),  (34)

where β₀ is the propagation factor for current propagating on thestructure. In the example of FIG. 3, I_(C) is the current that isdistributed along the vertical structure of the guided surface waveguideprobe 200 a.

For example, consider a feed network 209 that includes a low loss coil(e.g., a helical coil) at the bottom of the structure and a verticalfeed line conductor connected between the coil and the charge terminalT₁. The phase delay due to the coil (or helical delay line) isθ_(c)=β_(p)l_(C), with a physical length of l_(C) and a propagationfactor of

$\begin{matrix}{{\beta_{p} = {\frac{2\; \pi}{\lambda_{p}} = \frac{2\pi}{V_{f}\lambda_{0}}}},} & (35)\end{matrix}$

where V_(f) is the velocity factor on the structure, λ₀ is thewavelength at the supplied frequency, and λ_(p) is the propagationwavelength resulting from the velocity factor V_(f). The phase delay ismeasured relative to the ground (stake) current I₀.

In addition, the spatial phase delay along the length l_(w) of thevertical feed line conductor can be given by θ_(y)=β_(w)l_(w) whereβ_(w) is the propagation phase constant for the vertical feed lineconductor. In some implementations, the spatial phase delay may beapproximated by β_(y)=β_(w)h_(p), since the difference between thephysical height h_(p) of the guided surface waveguide probe 200 a andthe vertical feed line conductor length l_(w) is much less than awavelength at the supplied frequency (Δ₀). As a result, the total phasedelay through the coil and vertical feed line conductor isΦ=θ_(c)+θ_(y), and the current fed to the top of the coil from thebottom of the physical structure is

I _(C)(θ_(c)+θ_(y))=I ₀ e ^(jΦ),  (36)

with the total phase delay Φ measured relative to the ground (stake)current I₀. Consequently, the electrical effective height of a guidedsurface waveguide probe 200 can be approximated by

$\begin{matrix}{{h_{eff} = {{\frac{1}{I_{0}}{\int_{0}^{h_{p}}{I_{0}e^{j\; \Phi}{\cos \left( {\beta_{0}z} \right)}{dz}}}} \cong {h_{p}e^{j\; \Phi}}}},} & (37)\end{matrix}$

for the case where the physical height h_(p)<<Δ₀. The complex effectiveheight of a monopole, h_(eff)=h_(p) at an angle (or phase shift) of Φ,may be adjusted to cause the source fields to match a guided surfacewaveguide mode and cause a guided surface wave to be launched on thelossy conducting medium 203.

In the example of FIG. 5A, ray optics are used to illustrate the complexangle trigonometry of the incident electric field (E) having a complexBrewster angle of incidence (θ_(i,B)) at the Hankel crossover distance(R_(x)) 121. Recall from Equation (26) that, for a lossy conductingmedium, the Brewster angle is complex and specified by

$\begin{matrix}{{\tan \; \theta_{i,B}} = {\sqrt{ɛ_{r} - {j\frac{\sigma}{\omega \; ɛ_{o}}}} = {n.}}} & (38)\end{matrix}$

Electrically, the geometric parameters are related by the electricaleffective height (h_(eff)) of the charge terminal T₁ by

R _(x) tan ψ_(i,B) =R _(x) ×W=h _(eff) =h _(p) e ^(jΦ),  (39)

where ψ_(i,B)=(π/2)−θ_(i,B) is the Brewster angle measured from thesurface of the lossy conducting medium. To couple into the guidedsurface waveguide mode, the wave tilt of the electric field at theHankel crossover distance can be expressed as the ratio of theelectrical effective height and the Hankel crossover distance

$\begin{matrix}{\frac{h_{eff}}{R_{x}} = {{\tan \; \psi_{i,B}} = {W_{Rx}.}}} & (40)\end{matrix}$

Since both the physical height (h_(p)) and the Hankel crossover distance(R_(x)) are real quantities, the angle (Ψ) of the desired guided surfacewave tilt at the Hankel crossover distance (R_(x)) is equal to the phase(Φ) of the complex effective height (h_(eff)). This implies that byvarying the phase at the supply point of the coil, and thus the phaseshift in Equation (37), the phase, ϕ, of the complex effective heightcan be manipulated to match the angle of the wave tilt, ψ, of the guidedsurface waveguide mode at the Hankel crossover point 121: ϕ=ψ.

In FIG. 5A, a right triangle is depicted having an adjacent side oflength R_(x) along the lossy conducting medium surface and a complexBrewster angle ψ_(i,B) measured between a ray 124 extending between theHankel crossover point 121 at R_(x) and the center of the chargeterminal T₁, and the lossy conducting medium surface 127 between theHankel crossover point 121 and the charge terminal T₁. With the chargeterminal T₁ positioned at physical height h_(p) and excited with acharge having the appropriate phase delay Φ, the resulting electricfield is incident with the lossy conducting medium boundary interface atthe Hankel crossover distance R_(x), and at the Brewster angle. Underthese conditions, the guided surface waveguide mode can be excitedwithout reflection or substantially negligible reflection.

If the physical height of the charge terminal T₁ is decreased withoutchanging the phase shift Φ of the effective height (h_(eff)), theresulting electric field intersects the lossy conducting medium 203 atthe Brewster angle at a reduced distance from the guided surfacewaveguide probe 200. FIG. 6 graphically illustrates the effect ofdecreasing the physical height of the charge terminal T₁ on the distancewhere the electric field is incident at the Brewster angle. As theheight is decreased from h₃ through h₂ to h₁, the point where theelectric field intersects with the lossy conducting medium (e.g., theEarth) at the Brewster angle moves closer to the charge terminalposition. However, as Equation (39) indicates, the height H₁ (FIG. 3) ofthe charge terminal T₁ should be at or higher than the physical height(h_(p)) in order to excite the far-out component of the Hankel function.With the charge terminal T₁ positioned at or above the effective height(h_(eff)), the lossy conducting medium 203 can be illuminated at theBrewster angle of incidence (ψ_(i,B)=(π/2)−θ_(i,B)) at or beyond theHankel crossover distance (R_(x)) 121 as illustrated in FIG. 5A. Toreduce or minimize the bound charge on the charge terminal T₁, theheight should be at least four times the spherical diameter (orequivalent spherical diameter) of the charge terminal T₁ as mentionedabove.

A guided surface waveguide probe 200 can be configured to establish anelectric field having a wave tilt that corresponds to a waveilluminating the surface of the lossy conducting medium 203 at a complexBrewster angle, thereby exciting radial surface currents bysubstantially mode-matching to a guided surface wave mode at (or beyond)the Hankel crossover point 121 at R_(x).

Referring to FIG. 7, shown is a graphical representation of an exampleof a guided surface waveguide probe 200 b that includes a chargeterminal T₁. An AC source 212 acts as the excitation source for thecharge terminal T₁, which is coupled to the guided surface waveguideprobe 200 b through a feed network 209 (FIG. 3) comprising a coil 215such as, e.g., a helical coil. In other implementations, the AC source212 can be inductively coupled to the coil 215 through a primary coil.In some embodiments, an impedance matching network may be included toimprove and/or maximize coupling of the AC source 212 to the coil 215.

As shown in FIG. 7, the guided surface waveguide probe 200 b can includethe upper charge terminal T₁ (e.g., a sphere at height h_(p)) that ispositioned along a vertical axis z that is substantially normal to theplane presented by the lossy conducting medium 203. A second medium 206is located above the lossy conducting medium 203. The charge terminal T₁has a self-capacitance C_(T). During operation, charge Q₁ is imposed onthe terminal T₁ depending on the voltage applied to the terminal T₁ atany given instant.

In the example of FIG. 7, the coil 215 is coupled to a ground stake 218at a first end and to the charge terminal T₁ via a vertical feed lineconductor 221. In some implementations, the coil connection to thecharge terminal T₁ can be adjusted using a tap 224 of the coil 215 asshown in FIG. 7. The coil 215 can be energized at an operating frequencyby the AC source 212 through a tap 227 at a lower portion of the coil215. In other implementations, the AC source 212 can be inductivelycoupled to the coil 215 through a primary coil.

The construction and adjustment of the guided surface waveguide probe200 is based upon various operating conditions, such as the transmissionfrequency, conditions of the lossy conducting medium (e.g., soilconductivity σ and relative permittivity ε_(r)), and size of the chargeterminal T₁. The index of refraction can be calculated from Equations(10) and (11) as

n=√{square root over (ε_(r) −jx)},  (41)

where x=σ/ωε_(o) with ω=2πf. The conductivity σ and relativepermittivity ε_(r) can be determined through test measurements of thelossy conducting medium 203. The complex Brewster angle (θ_(i,B))measured from the surface normal can also be determined from Equation(26) as

θ_(i,B)=arctan(√{square root over (ε_(r) −jx)}),  (42)

or measured from the surface as shown in FIG. 5A as

$\begin{matrix}{\psi_{i,B} = {\frac{\pi}{2} - {\theta_{i,B}.}}} & (43)\end{matrix}$

The wave tilt at the Hankel crossover distance (W_(Rx)) can also befound using Equation (40).

The Hankel crossover distance can also be found by equating themagnitudes of Equations (20b) and (21) for −jγρ, and solving for R_(x)as illustrated by FIG. 4. The electrical effective height can then bedetermined from Equation (39) using the Hankel crossover distance andthe complex Brewster angle as

h _(eff) =h _(p) e ^(jΦ) =R _(x) tan ψ_(i,B).  (44)

As can be seen from Equation (44), the complex effective height(h_(eff)) includes a magnitude that is associated with the physicalheight (h_(p)) of the charge terminal T₁ and a phase delay (Φ) that isto be associated with the angle (Ψ) of the wave tilt at the Hankelcrossover distance (R_(x)). With these variables and the selected chargeterminal T₁ configuration, it is possible to determine the configurationof a guided surface waveguide probe 200.

With the charge terminal T₁ positioned at or above the physical height(h_(p)), the feed network 209 (FIG. 3) and/or the vertical feed lineconnecting the feed network to the charge terminal T₁ can be adjusted tomatch the phase (Φ) of the charge Q₁ on the charge terminal T₁ to theangle (Ψ) of the wave tilt (W). The size of the charge terminal T₁ canbe chosen to provide a sufficiently large surface for the charge Q₁imposed on the terminals. In general, it is desirable to make the chargeterminal T₁ as large as practical. The size of the charge terminal T₁should be large enough to avoid ionization of the surrounding air, whichcan result in electrical discharge or sparking around the chargeterminal.

The phase delay θ_(c), of a helically-wound coil can be determined fromMaxwell's equations as has been discussed by Corum, K. L. and J. F.Corum, “RF Coils, Helical Resonators and Voltage Magnification byCoherent Spatial Modes,” Microwave Review, Vol. 7, No. 2, September2001, pp. 36-45., which is incorporated herein by reference in itsentirety. For a helical coil with H/D>1, the ratio of the velocity ofpropagation (ν) of a wave along the coil's longitudinal axis to thespeed of light (c), or the “velocity factor,” is given by

$\begin{matrix}{{V_{f} = {\frac{\upsilon}{c} = \frac{1}{\sqrt{1 + {2\left( \frac{D}{s} \right)^{2.5}\left( \frac{D}{\lambda_{0}} \right)^{0.5}}}}}},} & (45)\end{matrix}$

where H is the axial length of the solenoidal helix, D is the coildiameter, N is the number of turns of the coil, s=H/N is theturn-to-turn spacing (or helix pitch) of the coil, and λ_(o) is thefree-space wavelength. Based upon this relationship, the electricallength, or phase delay, of the helical coil is given by

$\begin{matrix}{\theta_{c} = {{\beta_{p}H} = {{\frac{2\; \pi}{\lambda_{p}}H} = {\frac{2\; \pi}{V_{f}\lambda_{0}}{H.}}}}} & (46)\end{matrix}$

The principle is the same if the helix is wound spirally or is short andfat, but V_(f) and θ_(c) are easier to obtain by experimentalmeasurement. The expression for the characteristic (wave) impedance of ahelical transmission line has also been derived as

$\begin{matrix}{Z_{c} = {{\frac{60}{V_{f}}\left\lbrack {{\ln \left( \frac{V_{f}\lambda_{0}}{D} \right)} - 1.027} \right\rbrack}.}} & (47)\end{matrix}$

The spatial phase delay θ_(y) of the structure can be determined usingthe traveling wave phase delay of the vertical feed line conductor 221(FIG. 7). The capacitance of a cylindrical vertical conductor above aprefect ground plane can be expressed as

$\begin{matrix}{{C_{A} = {\frac{2\; \pi \; ɛ_{0}h_{w}}{{\ln \left( \frac{h}{a} \right)} - 1}{Farads}}},} & (48)\end{matrix}$

where h_(w) is the vertical length (or height) of the conductor and α isthe radius (in mks units). As with the helical coil, the traveling wavephase delay of the vertical feed line conductor can be given by

$\begin{matrix}{{\theta_{y} = {{\beta_{w}h_{w}} = {{\frac{2\; \pi}{\lambda_{w}}h_{w}} = {\frac{2\; \pi}{V_{w}\lambda_{0}}h_{w}}}}},} & (49)\end{matrix}$

where β_(w) is the propagation phase constant for the vertical feed lineconductor, h_(w) is the vertical length (or height) of the vertical feedline conductor, V_(w) is the velocity factor on the wire, λ₀ is thewavelength at the supplied frequency, and λ_(w) is the propagationwavelength resulting from the velocity factor V_(w). For a uniformcylindrical conductor, the velocity factor is a constant withV_(w)≈0.94, or in a range from about 0.93 to about 0.98. If the mast isconsidered to be a uniform transmission line, its average characteristicimpedance can be approximated by

$\begin{matrix}{{Z_{w} = {\frac{60}{V_{w}}\left\lbrack {{\ln \left( \frac{h_{w}}{a} \right)} - 1} \right\rbrack}},} & (50)\end{matrix}$

where V_(w)≈0.94 for a uniform cylindrical conductor and a is the radiusof the conductor. An alternative expression that has been employed inamateur radio literature for the characteristic impedance of asingle-wire feed line can be given by

$\begin{matrix}{Z_{w} = {138\mspace{14mu} {{\log \left( \frac{1.123\mspace{14mu} V_{w}\lambda_{0}}{2\; \pi \; a} \right)}.}}} & (51)\end{matrix}$

Equation (51) implies that Z_(w) for a single-wire feeder varies withfrequency. The phase delay can be determined based upon the capacitanceand characteristic impedance.

With a charge terminal T₁ positioned over the lossy conducting medium203 as shown in FIG. 3, the feed network 209 can be adjusted to excitethe charge terminal T₁ with the phase shift (Φ) of the complex effectiveheight (h_(eff)) equal to the angle (Ψ) of the wave tilt at the Hankelcrossover distance, or c, =T. When this condition is met, the electricfield produced by the charge oscillating Q₁ on the charge terminal T₁ iscoupled into a guided surface waveguide mode traveling along the surfaceof a lossy conducting medium 203. For example, if the Brewster angle(θ_(i,B)), the phase delay (θ_(y)) associated with the vertical feedline conductor 221 (FIG. 7), and the configuration of the coil 215 (FIG.7) are known, then the position of the tap 224 (FIG. 7) can bedetermined and adjusted to impose an oscillating charge Q₁ on the chargeterminal T₁ with phase Φ=Ψ. The position of the tap 224 may be adjustedto maximize coupling the traveling surface waves into the guided surfacewaveguide mode. Excess coil length beyond the position of the tap 224can be removed to reduce the capacitive effects. The vertical wireheight and/or the geometrical parameters of the helical coil may also bevaried.

The coupling to the guided surface waveguide mode on the surface of thelossy conducting medium 203 can be improved and/or optimized by tuningthe guided surface waveguide probe 200 for standing wave resonance withrespect to a complex image plane associated with the charge Q₁ on thecharge terminal T₁. By doing this, the performance of the guided surfacewaveguide probe 200 can be adjusted for increased and/or maximum voltage(and thus charge Q₁) on the charge terminal T₁. Referring back to FIG.3, the effect of the lossy conducting medium 203 in Region 1 can beexamined using image theory analysis.

Physically, an elevated charge Q₁ placed over a perfectly conductingplane attracts the free charge on the perfectly conducting plane, whichthen “piles up” in the region under the elevated charge Q₁. Theresulting distribution of “bound” electricity on the perfectlyconducting plane is similar to a bell-shaped curve. The superposition ofthe potential of the elevated charge Q₁, plus the potential of theinduced “piled up” charge beneath it, forces a zero equipotentialsurface for the perfectly conducting plane. The boundary value problemsolution that describes the fields in the region above the perfectlyconducting plane may be obtained using the classical notion of imagecharges, where the field from the elevated charge is superimposed withthe field from a corresponding “image” charge below the perfectlyconducting plane.

This analysis may also be used with respect to a lossy conducting medium203 by assuming the presence of an effective image charge Q₁′ beneaththe guided surface waveguide probe 200. The effective image charge Q₁′coincides with the charge Q₁ on the charge terminal T₁ about aconducting image ground plane 130, as illustrated in FIG. 3. However,the image charge Q₁′ is not merely located at some real depth and 180°out of phase with the primary source charge Q₁ on the charge terminalT₁, as they would be in the case of a perfect conductor. Rather, thelossy conducting medium 203 (e.g., a terrestrial medium) presents aphase shifted image. That is to say, the image charge Q₁′ is at acomplex depth below the surface (or physical boundary) of the lossyconducting medium 203. For a discussion of complex image depth,reference is made to Wait, J. R., “Complex Image Theory—Revisited,” IEEEAntennas and Propagation Maqazine, Vol. 33, No. 4, August 1991, pp.27-29, which is incorporated herein by reference in its entirety.

Instead of the image charge Q₁′ being at a depth that is equal to thephysical height (H₁) of the charge Q₁, the conducting image ground plane130 (representing a perfect conductor) is located at a complex depth ofz=−d/2 and the image charge Q₁′ appears at a complex depth (i.e., the“depth” has both magnitude and phase), given by −D₁=−(d/2+d/2+H₁)≠H₁.For vertically polarized sources over the Earth,

$\begin{matrix}{{d = {{\frac{2\sqrt{\gamma_{e}^{2} + k_{0}^{2}}}{\gamma_{e}^{2}} \approx \frac{2}{\gamma_{e}}} = {{d_{r} + {jd}_{i}} = {{d}{\angle\zeta}}}}},{where}} & (52) \\{{\gamma_{e}^{2} = {{j\; \omega \; \mu_{1}\sigma_{1}} - {\omega^{2}\mu_{1}ɛ_{1}}}},{and}} & (53) \\{{k_{0} = {\omega \sqrt{\mu_{0}ɛ_{0}}}},} & (54)\end{matrix}$

as indicated in Equation (12). The complex spacing of the image charge,in turn, implies that the external field will experience extra phaseshifts not encountered when the interface is either a dielectric or aperfect conductor. In the lossy conducting medium, the wave front normalis parallel to the tangent of the conducting image ground plane 130 atz=−d/2, and not at the boundary interface between Regions 1 and 2.

Consider the case illustrated in FIG. 8A where the lossy conductingmedium 203 is a finitely conducting Earth 133 with a physical boundary136. The finitely conducting Earth 133 may be replaced by a perfectlyconducting image ground plane 139 as shown in FIG. 8B, which is locatedat a complex depth z₁ below the physical boundary 136. This equivalentrepresentation exhibits the same impedance when looking down into theinterface at the physical boundary 136. The equivalent representation ofFIG. 8B can be modeled as an equivalent transmission line, as shown inFIG. 8C. The cross-section of the equivalent structure is represented asa (z-directed) end-loaded transmission line, with the impedance of theperfectly conducting image plane being a short circuit (z_(s)=0). Thedepth z₁ can be determined by equating the TEM wave impedance lookingdown at the Earth to an image ground plane impedance z_(in) seen lookinginto the transmission line of FIG. 8C.

In the case of FIG. 8A, the propagation constant and wave intrinsicimpedance in the upper region (air) 142 are

$\begin{matrix}{{\gamma_{o} = {{j\; \omega \sqrt{\mu_{o}ɛ_{o}}} = {0 + {j\; \beta_{o}}}}},{and}} & (55) \\{z_{o} = {\frac{j\; \omega \; \mu_{o}}{\gamma_{o}} = {\sqrt{\frac{\mu_{o}}{ɛ_{o}}}.}}} & (56)\end{matrix}$

In the lossy Earth 133, the propagation constant and wave intrinsicimpedance are

$\begin{matrix}{{\gamma_{e} = \sqrt{j\; \omega \; {\mu_{1}\left( {\sigma_{1} + {j\; \omega \; ɛ_{1}}} \right)}}},{and}} & (57) \\{Z_{e} = {\frac{j\; \omega \; \mu_{1}}{\gamma_{e}}.}} & (58)\end{matrix}$

For normal incidence, the equivalent representation of FIG. 8B isequivalent to a TEM transmission line whose characteristic impedance isthat of air (z_(o)), with propagation constant of γ_(o), and whoselength is z₁. As such, the image ground plane impedance Z_(in) seen atthe interface for the shorted transmission line of FIG. 8C is given by

Z _(in) =Z _(o) tan h(γ_(o) z ₁).  (59)

Equating the image ground plane impedance Z_(in) associated with theequivalent model of FIG. 8C to the normal incidence wave impedance ofFIG. 8A and solving for z₁ gives the distance to a short circuit (theperfectly conducting image ground plane 139) as

$\begin{matrix}{{z_{1} = {{\frac{1}{\gamma_{o}}{\tanh^{- 1}\left( \frac{Z_{e}}{Z_{o}} \right)}} = {{\frac{1}{\gamma_{o}}{\tanh^{- 1}\left( \frac{\gamma_{o}}{\gamma_{e}} \right)}} \approx \frac{1}{\gamma_{e}}}}},} & (60)\end{matrix}$

where only the first term of the series expansion for the inversehyperbolic tangent is considered for this approximation. Note that inthe air region 142, the propagation constant is γ_(o)=jβ_(o), soZ_(in)=jZ_(o) tan β_(o)z₁ (which is a purely imaginary quantity for areal z₁), but z_(e) is a complex value if σ≠0. Therefore, Z_(in)=Z_(e)only when z₁ is a complex distance.

Since the equivalent representation of FIG. 8B includes a perfectlyconducting image ground plane 139, the image depth for a charge orcurrent lying at the surface of the Earth (physical boundary 136) isequal to distance z₁ on the other side of the image ground plane 139, ord=2×z₁ beneath the Earth's surface (which is located at z=0). Thus, thedistance to the perfectly conducting image ground plane 139 can beapproximated by

$\begin{matrix}{d = {{2\; z_{1}} \approx {\frac{2}{\gamma_{e}}.}}} & (61)\end{matrix}$

Additionally, the “image charge” will be “equal and opposite” to thereal charge, so the potential of the perfectly conducting image groundplane 139 at depth z₁=−d/2 will be zero.

If a charge Q₁ is elevated a distance H₁ above the surface of the Earthas illustrated in FIG. 3, then the image charge Q₁′ resides at a complexdistance of D₁=d+H₁ below the surface, or a complex distance of d/2+H₁below the image ground plane 130. The guided surface waveguide probe 200b of FIG. 7 can be modeled as an equivalent single-wire transmissionline image plane model that can be based upon the perfectly conductingimage ground plane 139 of FIG. 8B. FIG. 9A shows an example of theequivalent single-wire transmission line image plane model, and FIG. 9Billustrates an example of the equivalent classic transmission linemodel, including the shorted transmission line of FIG. 8C.

In the equivalent image plane models of FIGS. 9A and 9B, Φ=θ_(y)+θ_(c)is the traveling wave phase delay of the guided surface waveguide probe200 referenced to Earth 133 (or the lossy conducting medium 203),θ_(c)=β_(p)H is the electrical length of the coil 215 (FIG. 7), ofphysical length H, expressed in degrees, θ_(y)=β_(w)h_(w) is theelectrical length of the vertical feed line conductor 221 (FIG. 7), ofphysical length h_(W), expressed in degrees, and θ_(d)=β_(o) d/2 is thephase shift between the image ground plane 139 and the physical boundary136 of the Earth 133 (or lossy conducting medium 203). In the example ofFIGS. 9A and 9B, Z_(W) is the characteristic impedance of the elevatedvertical feed line conductor 221 in ohms, Z_(c) is the characteristicimpedance of the coil 215 in ohms, and Z₀ is the characteristicimpedance of free space.

At the base of the guided surface waveguide probe 200, the impedanceseen “looking up” into the structure is Z_(↑)=Z_(base). With a loadimpedance of:

$\begin{matrix}{{Z_{L} = \frac{1}{j\; \omega \; C_{T}}},} & (62)\end{matrix}$

where C_(T) is the self-capacitance of the charge terminal T₁, theimpedance seen “looking up” into the vertical feed line conductor 221(FIG. 7) is given by:

$\begin{matrix}{{Z_{2} = {{Z_{W}\frac{Z_{L} + {Z_{w}\tan \; {h\left( {j\; \beta_{w}h_{w}} \right)}}}{Z_{w} + {Z_{L}\tan \; {h\left( {j\; \beta_{w}h_{w}} \right)}}}} = {Z_{W}\frac{Z_{L} + {Z_{w}\tan \; {h\left( {j\; \theta_{y}} \right)}}}{Z_{w} + {Z_{L}\tan \; {h\left( {j\; \theta_{y}} \right)}}}}}},} & (63)\end{matrix}$

and the impedance seen “looking up” into the coil 215 (FIG. 7) is givenby:

$\begin{matrix}{Z_{base} = {{Z_{c}\frac{Z_{2} + {Z_{c}\tan \; {h\left( {j\; \beta_{p}H} \right)}}}{Z_{c} + {Z_{2}\tan \; {h\left( {j\; \beta_{p}H} \right)}}}} = {Z_{c}{\frac{Z_{2} + {Z_{c}\tan \; {h\left( {j\; \theta_{c}} \right)}}}{Z_{c} + {Z_{2}\tan \; {h\left( {j\; \theta_{c}} \right)}}}.}}}} & (64)\end{matrix}$

At the base of the guided surface waveguide probe 200, the impedanceseen “looking down” into the lossy conducting medium 203 isZ_(↓)=Z_(in), which is given by:

$\begin{matrix}{{Z_{in} = {{Z_{o}\frac{Z_{s} + {Z_{o}\tan \; {h\left\lbrack {j\; {\beta_{o}\left( {d/2} \right)}} \right\rbrack}}}{Z_{o} + {Z_{s}\tan \; {h\left\lbrack {j\; {\beta_{o}\left( {d/2} \right)}} \right\rbrack}}}} = {Z_{o}\tan \; {h\left( {j\; \theta_{d}} \right)}}}},} & (65)\end{matrix}$

where Z_(s)=0.

Neglecting losses, the equivalent image plane model can be tuned toresonance when Z_(↓)+Z_(↑)=0 at the physical boundary 136. Or, in thelow loss case, X_(↓)+X_(↑)=0 at the physical boundary 136, where X isthe corresponding reactive component. Thus, the impedance at thephysical boundary 136 “looking up” into the guided surface waveguideprobe 200 is the conjugate of the impedance at the physical boundary 136“looking down” into the lossy conducting medium 203. By adjusting theload impedance Z_(L) of the charge terminal T₁ while maintaining thetraveling wave phase delay Φ equal to the angle of the media's wave tiltΨ, so that Φ=Ψ, which improves and/or maximizes coupling of the probe'selectric field to a guided surface waveguide mode along the surface ofthe lossy conducting medium 203 (e.g., Earth), the equivalent imageplane models of FIGS. 9A and 9B can be tuned to resonance with respectto the image ground plane 139. In this way, the impedance of theequivalent complex image plane model is purely resistive, whichmaintains a superposed standing wave on the probe structure thatmaximizes the voltage and elevated charge on terminal T₁, and byequations (1)-(3) and (16) maximizes the propagating surface wave.

It follows from the Hankel solutions, that the guided surface waveexcited by the guided surface waveguide probe 200 is an outwardpropagating traveling wave. The source distribution along the feednetwork 209 between the charge terminal T₁ and the ground stake 218 ofthe guided surface waveguide probe 200 (FIGS. 3 and 7) is actuallycomposed of a superposition of a traveling wave plus a standing wave onthe structure. With the charge terminal T₁ positioned at or above thephysical height h_(p), the phase delay of the traveling wave movingthrough the feed network 209 is matched to the angle of the wave tiltassociated with the lossy conducting medium 203. This mode-matchingallows the traveling wave to be launched along the lossy conductingmedium 203. Once the phase delay has been established for the travelingwave, the load impedance Z_(L) of the charge terminal T₁ is adjusted tobring the probe structure into standing wave resonance with respect tothe image ground plane (130 of FIG. 3 or 139 of FIG. 8), which is at acomplex depth of −d/2. In that case, the impedance seen from the imageground plane has zero reactance and the charge on the charge terminal T₁is maximized.

The distinction between the traveling wave phenomenon and standing wavephenomena is that (1) the phase delay of traveling waves (θ=βd) on asection of transmission line of length d (sometimes called a “delayline”) is due to propagation time delays; whereas (2) theposition-dependent phase of standing waves (which are composed offorward and backward propagating waves) depends on both the line lengthpropagation time delay and impedance transitions at interfaces betweenline sections of different characteristic impedances. In addition to thephase delay that arises due to the physical length of a section oftransmission line operating in sinusoidal steady-state, there is anextra reflection coefficient phase at impedance discontinuities that isdue to the ratio of Z_(oa)/Z_(ob), where Z_(oa) and Z_(ob) are thecharacteristic impedances of two sections of a transmission line suchas, e.g., a helical coil section of characteristic impedanceZ_(oa)=Z_(c) (FIG. 9B) and a straight section of vertical feed lineconductor of characteristic impedance Z_(ob)=Z_(w) (FIG. 9B).

As a result of this phenomenon, two relatively short transmission linesections of widely differing characteristic impedance may be used toprovide a very large phase shift. For example, a probe structurecomposed of two sections of transmission line, one of low impedance andone of high impedance, together totaling a physical length of, say,0.05λ, may be fabricated to provide a phase shift of 90° which isequivalent to a 0.25λ resonance. This is due to the large jump incharacteristic impedances. In this way, a physically short probestructure can be electrically longer than the two physical lengthscombined. This is illustrated in FIGS. 9A and 9B, where thediscontinuities in the impedance ratios provide large jumps in phase.The impedance discontinuity provides a substantial phase shift where thesections are joined together.

Referring to FIG. 10, shown is a flow chart 150 illustrating an exampleof adjusting a guided surface waveguide probe 200 (FIGS. 3 and 7) tosubstantially mode-match to a guided surface waveguide mode on thesurface of the lossy conducting medium, which launches a guided surfacetraveling wave along the surface of a lossy conducting medium 203 (FIG.3). Beginning with 153, the charge terminal T₁ of the guided surfacewaveguide probe 200 is positioned at a defined height above a lossyconducting medium 203. Utilizing the characteristics of the lossyconducting medium 203 and the operating frequency of the guided surfacewaveguide probe 200, the Hankel crossover distance can also be found byequating the magnitudes of Equations (20b) and (21) for −jγρ, andsolving for R_(x) as illustrated by FIG. 4. The complex index ofrefraction (n) can be determined using Equation (41), and the complexBrewster angle (θ_(i,B)) can then be determined from Equation (42). Thephysical height (h_(p)) of the charge terminal T₁ can then be determinedfrom Equation (44). The charge terminal T₁ should be at or higher thanthe physical height (h_(p)) in order to excite the far-out component ofthe Hankel function. This height relationship is initially consideredwhen launching surface waves. To reduce or minimize the bound charge onthe charge terminal T₁, the height should be at least four times thespherical diameter (or equivalent spherical diameter) of the chargeterminal T₁.

At 156, the electrical phase delay Φ of the elevated charge Q₁ on thecharge terminal T₁ is matched to the complex wave tilt angle Ψ. Thephase delay (θ_(c)) of the helical coil and/or the phase delay (θ_(y))of the vertical feed line conductor can be adjusted to make Φ equal tothe angle (Ψ) of the wave tilt (W). Based on Equation (31), the angle(Ψ) of the wave tilt can be determined from:

$\begin{matrix}{W = {\frac{E_{\rho}}{E_{z}} = {\frac{1}{\tan \; \theta_{i,B}} = {\frac{1}{n} = {{W}{e^{j\; \Psi}.}}}}}} & (66)\end{matrix}$

The electrical phase Φ can then be matched to the angle of the wavetilt. This angular (or phase) relationship is next considered whenlaunching surface waves. For example, the electrical phase delayΦ=θ_(c)+θ_(y) can be adjusted by varying the geometrical parameters ofthe coil 215 (FIG. 7) and/or the length (or height) of the vertical feedline conductor 221 (FIG. 7). By matching Φ=Ψ, an electric field can beestablished at or beyond the Hankel crossover distance (R_(x)) with acomplex Brewster angle at the boundary interface to excite the surfacewaveguide mode and launch a traveling wave along the lossy conductingmedium 203.

Next at 159, the load impedance of the charge terminal T₁ is tuned toresonate the equivalent image plane model of the guided surfacewaveguide probe 200. The depth (d/2) of the conducting image groundplane 139 of FIGS. 9A and 9B (or 130 of FIG. 3) can be determined usingEquations (52), (53) and (54) and the values of the lossy conductingmedium 203 (e.g., the Earth), which can be measured. Using that depth,the phase shift (θ_(d)) between the image ground plane 139 and thephysical boundary 136 of the lossy conducting medium 203 can bedetermined using θ_(d)=β_(o) d/2. The impedance (Z_(in)) as seen“looking down” into the lossy conducting medium 203 can then bedetermined using Equation (65). This resonance relationship can beconsidered to maximize the launched surface waves.

Based upon the adjusted parameters of the coil 215 and the length of thevertical feed line conductor 221, the velocity factor, phase delay, andimpedance of the coil 215 and vertical feed line conductor 221 can bedetermined using Equations (45) through (51). In addition, theself-capacitance (C_(T)) of the charge terminal T₁ can be determinedusing, e.g., Equation (24). The propagation factor (β_(p)) of the coil215 can be determined using Equation (35) and the propagation phaseconstant (β_(w)) for the vertical feed line conductor 221 can bedetermined using Equation (49). Using the self-capacitance and thedetermined values of the coil 215 and vertical feed line conductor 221,the impedance (Z_(base)) of the guided surface waveguide probe 200 asseen “looking up” into the coil 215 can be determined using Equations(62), (63) and (64).

The equivalent image plane model of the guided surface waveguide probe200 can be tuned to resonance by adjusting the load impedance Z_(L) suchthat the reactance component X_(base) of Z_(base) cancels out thereactance component X_(in) of Z_(in), or X_(base)+X_(in)=0. Thus, theimpedance at the physical boundary 136 “looking up” into the guidedsurface waveguide probe 200 is the conjugate of the impedance at thephysical boundary 136 “looking down” into the lossy conducting medium203. The load impedance Z_(L) can be adjusted by varying the capacitance(C_(T)) of the charge terminal T₁ without changing the electrical phasedelay Φ=θ_(c)+θ_(y) of the charge terminal T₁. An iterative approach maybe taken to tune the load impedance Z_(L) for resonance of theequivalent image plane model with respect to the conducting image groundplane 139 (or 130). In this way, the coupling of the electric field to aguided surface waveguide mode along the surface of the lossy conductingmedium 203 (e.g., Earth) can be improved and/or maximized.

This may be better understood by illustrating the situation with anumerical example. Consider a guided surface waveguide probe 200comprising a top-loaded vertical stub of physical height h_(p) with acharge terminal T₁ at the top, where the charge terminal T₁ is excitedthrough a helical coil and vertical feed line conductor at anoperational frequency (f_(o)) of 1.85 MHz. With a height (H₁) of 16 feetand the lossy conducting medium 203 (e.g., Earth) having a relativepermittivity of ε_(r)=15 and a conductivity of σ₁=0.010 mhos/m, severalsurface wave propagation parameters can be calculated for f_(o)=1.850MHz. Under these conditions, the Hankel crossover distance can be foundto be R_(x)=54.5 feet with a physical height of h_(p)=5.5 feet, which iswell below the actual height of the charge terminal T₁. While a chargeterminal height of H₁=5.5 feet could have been used, the taller probestructure reduced the bound capacitance, permitting a greater percentageof free charge on the charge terminal T₁ providing greater fieldstrength and excitation of the traveling wave.

The wave length can be determined as:

$\begin{matrix}{{\lambda_{o} = {\frac{c}{f_{o}} = {162.162\mspace{14mu} {meters}}}},} & (67)\end{matrix}$

where c is the speed of light. The complex index of refraction is:

n=√{square root over (ε_(r) −jx)}=7.529−j6.546,  (68)

from Equation (41), where x=σ₁/ωε_(o) with ω=2πf_(o), and the complexBrewster angle is:

θ_(i,B)=arctan(√{square root over (ε_(r) −jx)})=85.6−j3.744°.  (69)

from Equation (42). Using Equation (66), the wave tilt values can bedetermined to be:

$\begin{matrix}{W = {\frac{1}{\tan \; \theta_{i,B}} = {\frac{1}{n} = {{{W}e^{j\; \Psi}} = {0.101\; {e^{j\; {40.614 \circ}}.}}}}}} & (70)\end{matrix}$

Thus, the helical coil can be adjusted to match Φ=Ψ=40.614°

The velocity factor of the vertical feed line conductor (approximated asa uniform cylindrical conductor with a diameter of 0.27 inches) can begiven as V_(w)≈0.93. Since h_(p)<<λ_(o), the propagation phase constantfor the vertical feed line conductor can be approximated as:

$\begin{matrix}{\beta_{w} = {\frac{2\; \pi}{\lambda_{w}} = {\frac{2\; \pi}{V_{w}\lambda_{0}} = {0.042\mspace{14mu} {m^{- 1}.}}}}} & (71)\end{matrix}$

From Equation (49) the phase delay of the vertical feed line conductoris:

θ_(y)=β_(w) h _(w)≈β_(w) h _(p)=11.640°.  (72)

By adjusting the phase delay of the helical coil so thatθ_(c)=28.974°=40.614−11.640°, Φ will equal Ψ to match the guided surfacewaveguide mode. To illustrate the relationship between Φ and Ψ, FIG. 11shows a plot of both over a range of frequencies. As both Φ and Ψ arefrequency dependent, it can be seen that their respective curves crossover each other at approximately 1.85 MHz.

For a helical coil having a conductor diameter of 0.0881 inches, a coildiameter (D) of 30 inches and a turn-to-turn spacing (s) of 4 inches,the velocity factor for the coil can be determined using Equation (45)as:

$\begin{matrix}{{V_{f} = {\frac{1}{\sqrt{1 + {20\left( \frac{D}{s} \right)^{2.5}\left( \frac{D}{\lambda_{o}} \right)^{0.5}}}} = 0.069}},} & (73)\end{matrix}$

and the propagation factor from Equation (35) is:

$\begin{matrix}{\beta_{p} = {\frac{2\; \pi}{V_{f}\lambda_{0}} = {0.564\mspace{14mu} {m^{- 1}.}}}} & (74)\end{matrix}$

With θ_(c)=28.974°, the axial length of the solenoidal helix (H) can bedetermined using Equation (46) such that:

$\begin{matrix}{H = {\frac{\theta_{c}}{\beta_{p}} = {35.2732\mspace{14mu} {{inches}.}}}} & (75)\end{matrix}$

This height determines the location on the helical coil where thevertical feed line conductor is connected, resulting in a coil with8.818 turns (N=H/s).

With the traveling wave phase delay of the coil and vertical feed lineconductor adjusted to match the wave tilt angle (Φ=θ_(c)+θ_(y)=Ψ), theload impedance (Z_(L)) of the charge terminal T₁ can be adjusted forstanding wave resonance of the equivalent image plane model of theguided surface wave probe 200. From the measured permittivity,conductivity and permeability of the Earth, the radial propagationconstant can be determined using Equation (57)

γ_(e)=√{square root over (jωu ₁(σ₁ +jωε ₁))}=0.25+j0.292m ⁻¹,  (76)

And the complex depth of the conducting image ground plane can beapproximated from Equation (52) as:

$\begin{matrix}{{{d \approx \frac{2}{\gamma_{e}}} = {3.364 + {j\mspace{14mu} 3.963\mspace{14mu} {meters}}}},} & (77)\end{matrix}$

with a corresponding phase shift between the conducting image groundplane and the physical boundary of the Earth given by:

θ_(d)=β_(o)(d/2)=4.015−j4.73°.  (78)

Using Equation (65), the impedance seen “looking down” into the lossyconducting medium 203 (i.e., Earth) can be determined as:

Z _(in) =Z _(o) tan h(jθ _(d))=R _(in) +jX _(in)=31.191+j26.27ohms.  (79)

By matching the reactive component (X_(in)) seen “looking down” into thelossy conducting medium 203 with the reactive component (X_(base)) seen“looking up” into the guided surface wave probe 200, the coupling intothe guided surface waveguide mode may be maximized. This can beaccomplished by adjusting the capacitance of the charge terminal T₁without changing the traveling wave phase delays of the coil andvertical feed line conductor. For example, by adjusting the chargeterminal capacitance (C_(T)) to 61.8126 pF, the load impedance fromEquation (62) is:

$\begin{matrix}{{Z_{L} = {\frac{1}{j\; \omega \; C_{T}} = {{- j}\mspace{14mu} 1392\mspace{14mu} {ohms}}}},} & (80)\end{matrix}$

and the reactive components at the boundary are matched.

Using Equation (51), the impedance of the vertical feed line conductor(having a diameter (2a) of 0.27 inches) is given as

$\begin{matrix}{{Z_{w} = {{138\mspace{14mu} {\log \left( \frac{1.123\mspace{14mu} V_{w}\lambda_{0}}{2\; \pi \; a} \right)}} = {537.534\mspace{14mu} {ohms}}}},} & (81)\end{matrix}$

and the impedance seen “looking up” into the vertical feed lineconductor is given by Equation (63) as:

$\begin{matrix}{Z_{2} = {{Z_{W}\frac{Z_{L} + {Z_{w}{\tanh \left( {j\; \theta_{y}} \right)}}}{Z_{w} + {Z_{L}\tanh \; \left( {j\; \theta_{y}} \right)}}} = {{- j}\mspace{14mu} 835.438\mspace{14mu} {{ohms}.}}}} & (82)\end{matrix}$

Using Equation (47), the characteristic impedance of the helical coil isgiven as

$\begin{matrix}{{Z_{c} = {{\frac{60}{V_{f}}\left\lbrack {{\ln\left( \frac{V_{f}\lambda_{0}}{D} \right)} - 1.027} \right\rbrack} = {1446\mspace{14mu} {ohms}}}},} & (83)\end{matrix}$

and the impedance seen “looking up” into the coil at the base is givenby Equation (64) as:

$\begin{matrix}{Z_{base} = {{Z_{c}\frac{Z_{2} + {Z_{c}{\tanh \left( {j\; \theta_{c}} \right)}}}{Z_{c} + {Z_{2}{\tanh \left( {j\; \theta_{c}} \right)}}}} = {{- j}\mspace{14mu} 26.271\mspace{14mu} {{ohms}.}}}} & (84)\end{matrix}$

When compared to the solution of Equation (79), it can be seen that thereactive components are opposite and approximately equal, and thus areconjugates of each other. Thus, the impedance (Z_(ip)) seen “looking up”into the equivalent image plane model of FIGS. 9A and 9B from theperfectly conducting image ground plane is only resistive orZ_(ip)=R+j0.

When the electric fields produced by a guided surface waveguide probe200 (FIG. 3) are established by matching the traveling wave phase delayof the feed network to the wave tilt angle and the probe structure isresonated with respect to the perfectly conducting image ground plane atcomplex depth z=−d/2, the fields are substantially mode-matched to aguided surface waveguide mode on the surface of the lossy conductingmedium, a guided surface traveling wave is launched along the surface ofthe lossy conducting medium. As illustrated in FIG. 1, the guided fieldstrength curve 103 of the guided electromagnetic field has acharacteristic exponential decay of e^(−ad)/√{square root over (d)} andexhibits a distinctive knee 109 on the log-log scale.

In summary, both analytically and experimentally, the traveling wavecomponent on the structure of the guided surface waveguide probe 200 hasa phase delay (Φ) at its upper terminal that matches the angle (Ψ) ofthe wave tilt of the surface traveling wave (Φ=Ψ). Under this condition,the surface waveguide may be considered to be “mode-matched”.Furthermore, the resonant standing wave component on the structure ofthe guided surface waveguide probe 200 has a V_(MAX) at the chargeterminal T₁ and a V_(MIN) down at the image plane 139 (FIG. 8B) whereZ_(ip)=R_(ip)+j0 at a complex depth of z=−d/2, not at the connection atthe physical boundary 136 of the lossy conducting medium 203 (FIG. 8B).Lastly, the charge terminal T₁ is of sufficient height H₁ of FIG. 3(h≥R_(x) tan ψ_(i,B)) so that electromagnetic waves incident onto thelossy conducting medium 203 at the complex Brewster angle do so out at adistance (≥R_(x)) where the 1/√{square root over (r)} term ispredominant. Receive circuits can be utilized with one or more guidedsurface waveguide probes to facilitate wireless transmission and/orpower delivery systems.

Referring back to FIG. 3, operation of a guided surface waveguide probe200 may be controlled to adjust for variations in operational conditionsassociated with the guided surface waveguide probe 200. For example, anadaptive probe control system 230 can be used to control the feednetwork 209 and/or the charge terminal T₁ to control the operation ofthe guided surface waveguide probe 200. Operational conditions caninclude, but are not limited to, variations in the characteristics ofthe lossy conducting medium 203 (e.g., conductivity σ and relativepermittivity ε_(r)), variations in field strength and/or variations inloading of the guided surface waveguide probe 200. As can be seen fromEquations (31), (41) and (42), the index of refraction (n), the complexBrewster angle (θ_(i,B)), and the wave tilt (|W|e^(jΨ)) can be affectedby changes in soil conductivity and permittivity resulting from, e.g.,weather conditions.

Equipment such as, e.g., conductivity measurement probes, permittivitysensors, ground parameter meters, field meters, current monitors and/orload receivers can be used to monitor for changes in the operationalconditions and provide information about current operational conditionsto the adaptive probe control system 230. The probe control system 230can then make one or more adjustments to the guided surface waveguideprobe 200 to maintain specified operational conditions for the guidedsurface waveguide probe 200. For instance, as the moisture andtemperature vary, the conductivity of the soil will also vary.Conductivity measurement probes and/or permittivity sensors may belocated at multiple locations around the guided surface waveguide probe200. Generally, it would be desirable to monitor the conductivity and/orpermittivity at or about the Hankel crossover distance R_(x) for theoperational frequency. Conductivity measurement probes and/orpermittivity sensors may be located at multiple locations (e.g., in eachquadrant) around the guided surface waveguide probe 200.

The conductivity measurement probes and/or permittivity sensors can beconfigured to evaluate the conductivity and/or permittivity on aperiodic basis and communicate the information to the probe controlsystem 230. The information may be communicated to the probe controlsystem 230 through a network such as, but not limited to, a LAN, WLAN,cellular network, or other appropriate wired or wireless communicationnetwork. Based upon the monitored conductivity and/or permittivity, theprobe control system 230 may evaluate the variation in the index ofrefraction (n), the complex Brewster angle (θ_(i,B)), and/or the wavetilt (|W|e^(jΨ)) and adjust the guided surface waveguide probe 200 tomaintain the phase delay (Φ) of the feed network 209 equal to the wavetilt angle (Ψ) and/or maintain resonance of the equivalent image planemodel of the guided surface waveguide probe 200. This can beaccomplished by adjusting, e.g., θ_(y), θ_(c) and/or C_(T). Forinstance, the probe control system 230 can adjust the self-capacitanceof the charge terminal T₁ and/or the phase delay (θy, θ_(c)) applied tothe charge terminal T₁ to maintain the electrical launching efficiencyof the guided surface wave at or near its maximum. For example, theself-capacitance of the charge terminal T₁ can be varied by changing thesize of the terminal. The charge distribution can also be improved byincreasing the size of the charge terminal T₁, which can reduce thechance of an electrical discharge from the charge terminal T₁. In otherembodiments, the charge terminal T₁ can include a variable inductancethat can be adjusted to change the load impedance Z_(L). The phaseapplied to the charge terminal T₁ can be adjusted by varying the tapposition on the coil 215 (FIG. 7), and/or by including a plurality ofpredefined taps along the coil 215 and switching between the differentpredefined tap locations to maximize the launching efficiency.

Field or field strength (FS) meters may also be distributed about theguided surface waveguide probe 200 to measure field strength of fieldsassociated with the guided surface wave. The field or FS meters can beconfigured to detect the field strength and/or changes in the fieldstrength (e.g., electric field strength) and communicate thatinformation to the probe control system 230. The information may becommunicated to the probe control system 230 through a network such as,but not limited to, a LAN, WLAN, cellular network, or other appropriatecommunication network. As the load and/or environmental conditionschange or vary during operation, the guided surface waveguide probe 200may be adjusted to maintain specified field strength(s) at the FS meterlocations to ensure appropriate power transmission to the receivers andthe loads they supply.

For example, the phase delay (Ψ=θ_(y)+θ_(c)) applied to the chargeterminal T₁ can be adjusted to match the wave tilt angle (Ψ). Byadjusting one or both phase delays, the guided surface waveguide probe200 can be adjusted to ensure the wave tilt corresponds to the complexBrewster angle. This can be accomplished by adjusting a tap position onthe coil 215 (FIG. 7) to change the phase delay supplied to the chargeterminal T₁. The voltage level supplied to the charge terminal T₁ canalso be increased or decreased to adjust the electric field strength.This may be accomplished by adjusting the output voltage of theexcitation source 212 or by adjusting or reconfiguring the feed network209. For instance, the position of the tap 227 (FIG. 7) for the ACsource 212 can be adjusted to increase the voltage seen by the chargeterminal T₁. Maintaining field strength levels within predefined rangescan improve coupling by the receivers, reduce ground current losses, andavoid interference with transmissions from other guided surfacewaveguide probes 200.

The probe control system 230 can be implemented with hardware, firmware,software executed by hardware, or a combination thereof. For example,the probe control system 230 can include processing circuitry includinga processor and a memory, both of which can be coupled to a localinterface such as, for example, a data bus with an accompanyingcontrol/address bus as can be appreciated by those with ordinary skillin the art. A probe control application may be executed by the processorto adjust the operation of the guided surface waveguide probe 200 basedupon monitored conditions. The probe control system 230 can also includeone or more network interfaces for communicating with the variousmonitoring devices. Communications can be through a network such as, butnot limited to, a LAN, WLAN, cellular network, or other appropriatecommunication network. The probe control system 230 may comprise, forexample, a computer system such as a server, desktop computer, laptop,or other system with like capability.

Referring back to the example of FIG. 5A, the complex angle trigonometryis shown for the ray optic interpretation of the incident electric field(E) of the charge terminal T₁ with a complex Brewster angle (θ_(i,B)) atthe Hankel crossover distance (R_(x)). Recall that, for a lossyconducting medium, the Brewster angle is complex and specified byequation (38). Electrically, the geometric parameters are related by theelectrical effective height (h_(eff)) of the charge terminal T₁ byequation (39). Since both the physical height (h_(p)) and the Hankelcrossover distance (R_(x)) are real quantities, the angle of the desiredguided surface wave tilt at the Hankel crossover distance (W_(Rx)) isequal to the phase (Ψ) of the complex effective height (h_(eff)). Withthe charge terminal T₁ positioned at the physical height h_(p) andexcited with a charge having the appropriate phase Ψ, the resultingelectric field is incident with the lossy conducting medium boundaryinterface at the Hankel crossover distance R_(x), and at the Brewsterangle. Under these conditions, the guided surface waveguide mode can beexcited without reflection or substantially negligible reflection.

However, Equation (39) means that the physical height of the guidedsurface waveguide probe 200 can be relatively small. While this willexcite the guided surface waveguide mode, this can result in an undulylarge bound charge with little free charge. To compensate, the chargeterminal T₁ can be raised to an appropriate elevation to increase theamount of free charge. As one example rule of thumb, the charge terminalT₁ can be positioned at an elevation of about 4-5 times (or more) theeffective diameter of the charge terminal T₁. FIG. 6 illustrates theeffect of raising the charge terminal T₁ above the physical height(h_(p)) shown in FIG. 5A. The increased elevation causes the distance atwhich the wave tilt is incident with the lossy conductive medium to movebeyond the Hankel crossover point 121 (FIG. 5A). To improve coupling inthe guided surface waveguide mode, and thus provide for a greaterlaunching efficiency of the guided surface wave, a lower compensationterminal T₂ can be used to adjust the total effective height (h_(TE)) ofthe charge terminal T₁ such that the wave tilt at the Hankel crossoverdistance is at the Brewster angle.

Referring to FIG. 12, shown is an example of a guided surface waveguideprobe 200 c that includes an elevated charge terminal T₁ and a lowercompensation terminal T₂ that are arranged along a vertical axis z thatis normal to a plane presented by the lossy conducting medium 203. Inthis respect, the charge terminal T₁ is placed directly above thecompensation terminal T₂ although it is possible that some otherarrangement of two or more charge and/or compensation terminals TN canbe used. The guided surface waveguide probe 200 c is disposed above alossy conducting medium 203 according to an embodiment of the presentdisclosure. The lossy conducting medium 203 makes up Region 1 with asecond medium 206 that makes up Region 2 sharing a boundary interfacewith the lossy conducting medium 203.

The guided surface waveguide probe 200 c includes a feed network 209that couples an excitation source 212 to the charge terminal T₁ and thecompensation terminal T₂. According to various embodiments, charges Q₁and Q₂ can be imposed on the respective charge and compensationterminals T₁ and T₂, depending on the voltages applied to terminals T₁and T₂ at any given instant. I₁ is the conduction current feeding thecharge Q₁ on the charge terminal T₁ via the terminal lead, and I₂ is theconduction current feeding the charge Q₂ on the compensation terminal T₂via the terminal lead.

According to the embodiment of FIG. 12, the charge terminal T₁ ispositioned over the lossy conducting medium 203 at a physical height H₁,and the compensation terminal T₂ is positioned directly below T₁ alongthe vertical axis z at a physical height H₂, where H₂ is less than H₁.The height h of the transmission structure may be calculated as h=H₁−H₂.The charge terminal T₁ has an isolated (or self) capacitance C₁, and thecompensation terminal T₂ has an isolated (or self) capacitance C₂. Amutual capacitance CM can also exist between the terminals T₁ and T₂depending on the distance therebetween. During operation, charges Q₁ andQ₂ are imposed on the charge terminal T₁ and the compensation terminalT₂, respectively, depending on the voltages applied to the chargeterminal T₁ and the compensation terminal T₂ at any given instant.

Referring next to FIG. 13, shown is a ray optics interpretation of theeffects produced by the elevated charge Q₁ on charge terminal T₁ andcompensation terminal T₂ of FIG. 12. With the charge terminal T₁elevated to a height where the ray intersects with the lossy conductivemedium at the Brewster angle at a distance greater than the Hankelcrossover point 121 as illustrated by line 163, the compensationterminal T₂ can be used to adjust h_(TE) by compensating for theincreased height. The effect of the compensation terminal T₂ is toreduce the electrical effective height of the guided surface waveguideprobe (or effectively raise the lossy medium interface) such that thewave tilt at the Hankel crossover distance is at the Brewster angle asillustrated by line 166.

The total effective height can be written as the superposition of anupper effective height (h_(UE)) associated with the charge terminal T₁and a lower effective height (h_(LE)) associated with the compensationterminal T₂ such that

h _(TE) =h _(UE) +h _(LE) =h _(p) e ^(j(βh) ^(p) ^(+Φ) ^(U) ⁾ +h _(d) e^(j(βh) ^(d) ^(+Φ) ^(L) ⁾ =R _(x) ×W,  (85)

where Ψ_(U) is the phase delay applied to the upper charge terminal T₁,Ψ_(L) is the phase delay applied to the lower compensation terminal T₂,β=2π/λ_(p) is the propagation factor from Equation (35), h_(p) is thephysical height of the charge terminal T₁ and h_(d) is the physicalheight of the compensation terminal T₂. If extra lead lengths are takeninto consideration, they can be accounted for by adding the chargeterminal lead length z to the physical height h_(p) of the chargeterminal T₁ and the compensation terminal lead length y to the physicalheight h_(d) of the compensation terminal T₂ as shown in

h _(TE)=(h _(p) +Z)e ^(j(β(h) ^(p) ^(+z)+Φ) ^(U) ⁾+(h _(d) +y)e ^(j(β(h)^(d) ^(+y)+Φ) ^(L) ⁾ =R _(x) ×W.  (86)

The lower effective height can be used to adjust the total effectiveheight (h_(TE)) to equal the complex effective height (h_(eff)) of FIG.5A.

Equations (85) or (86) can be used to determine the physical height ofthe lower disk of the compensation terminal T₂ and the phase angles tofeed the terminals in order to obtain the desired wave tilt at theHankel crossover distance. For example, Equation (86) can be rewrittenas the phase shift applied to the charge terminal T₁ as a function ofthe compensation terminal height (h_(d)) to give

$\begin{matrix}{{\Phi_{U}\left( h_{d} \right)} = {{- {\beta \left( {h_{p} + z} \right)}} - {j\mspace{14mu} {{\ln\left( \frac{{R_{x} \times W} - {\left( {h_{d} + y} \right)e^{j{({{\beta \; h_{d}} + {\beta \; y} + \Phi_{L}})}}}}{\left( {h_{p} + z} \right)} \right)}.}}}} & (87)\end{matrix}$

To determine the positioning of the compensation terminal T₂, therelationships discussed above can be utilized. First, the totaleffective height (h_(TE)) is the superposition of the complex effectiveheight (h_(UE)) of the upper charge terminal T₁ and the complexeffective height (h_(LE)) of the lower compensation terminal T₂ asexpressed in Equation (86). Next, the tangent of the angle of incidencecan be expressed geometrically as

$\begin{matrix}{{{\tan \; \psi_{E}} = \frac{h_{TE}}{R_{x}}},} & (88)\end{matrix}$

which is equal to the definition of the wave tilt, W. Finally, given thedesired Hankel crossover distance R_(x), the h_(TE) can be adjusted tomake the wave tilt of the incident ray match the complex Brewster angleat the Hankel crossover point 121. This can be accomplished by adjustingh_(p), Ψ_(U), and/or h_(d).

These concepts may be better understood when discussed in the context ofan example of a guided surface waveguide probe. Referring to FIG. 14,shown is a graphical representation of an example of a guided surfacewaveguide probe 200 d including an upper charge terminal T₁ (e.g., asphere at height h_(T)) and a lower compensation terminal T₂ (e.g., adisk at height h_(d)) that are positioned along a vertical axis z thatis substantially normal to the plane presented by the lossy conductingmedium 203. During operation, charges Q₁ and Q₂ are imposed on thecharge and compensation terminals T₁ and T₂, respectively, depending onthe voltages applied to the terminals T₁ and T₂ at any given instant.

An AC source 212 acts as the excitation source for the charge terminalT₁, which is coupled to the guided surface waveguide probe 200 d througha feed network 209 comprising a coil 215 such as, e.g., a helical coil.The AC source 212 can be connected across a lower portion of the coil215 through a tap 227, as shown in FIG. 14, or can be inductivelycoupled to the coil 215 by way of a primary coil. The coil 215 can becoupled to a ground stake 218 at a first end and the charge terminal T₁at a second end. In some implementations, the connection to the chargeterminal T₁ can be adjusted using a tap 224 at the second end of thecoil 215. The compensation terminal T₂ is positioned above andsubstantially parallel with the lossy conducting medium 203 (e.g., theground or Earth), and energized through a tap 233 coupled to the coil215. An ammeter 236 located between the coil 215 and ground stake 218can be used to provide an indication of the magnitude of the currentflow (I₀) at the base of the guided surface waveguide probe.Alternatively, a current clamp may be used around the conductor coupledto the ground stake 218 to obtain an indication of the magnitude of thecurrent flow (I₀).

In the example of FIG. 14, the coil 215 is coupled to a ground stake 218at a first end and the charge terminal T₁ at a second end via a verticalfeed line conductor 221. In some implementations, the connection to thecharge terminal T₁ can be adjusted using a tap 224 at the second end ofthe coil 215 as shown in FIG. 14. The coil 215 can be energized at anoperating frequency by the AC source 212 through a tap 227 at a lowerportion of the coil 215. In other implementations, the AC source 212 canbe inductively coupled to the coil 215 through a primary coil. Thecompensation terminal T₂ is energized through a tap 233 coupled to thecoil 215. An ammeter 236 located between the coil 215 and ground stake218 can be used to provide an indication of the magnitude of the currentflow at the base of the guided surface waveguide probe 200 d.Alternatively, a current clamp may be used around the conductor coupledto the ground stake 218 to obtain an indication of the magnitude of thecurrent flow. The compensation terminal T₂ is positioned above andsubstantially parallel with the lossy conducting medium 203 (e.g., theground).

In the example of FIG. 14, the connection to the charge terminal T₁located on the coil 215 above the connection point of tap 233 for thecompensation terminal T₂. Such an adjustment allows an increased voltage(and thus a higher charge Q₁) to be applied to the upper charge terminalT₁. In other embodiments, the connection points for the charge terminalT₁ and the compensation terminal T₂ can be reversed. It is possible toadjust the total effective height (h_(TE)) of the guided surfacewaveguide probe 200 d to excite an electric field having a guidedsurface wave tilt at the Hankel crossover distance R_(x). The Hankelcrossover distance can also be found by equating the magnitudes ofequations (20b) and (21) for −jγρ, and solving for R_(x) as illustratedby FIG. 4. The index of refraction (n), the complex Brewster angle(θ_(i,B) and ψ_(i,B)), the wave tilt (|W|e^(jΨ)) and the complexeffective height (h_(eff)=h_(p)e^(jΦ)) can be determined as describedwith respect to Equations (41)-(44) above.

With the selected charge terminal T₁ configuration, a spherical diameter(or the effective spherical diameter) can be determined. For example, ifthe charge terminal T₁ is not configured as a sphere, then the terminalconfiguration may be modeled as a spherical capacitance having aneffective spherical diameter. The size of the charge terminal T₁ can bechosen to provide a sufficiently large surface for the charge Q₁ imposedon the terminals. In general, it is desirable to make the chargeterminal T₁ as large as practical. The size of the charge terminal T₁should be large enough to avoid ionization of the surrounding air, whichcan result in electrical discharge or sparking around the chargeterminal. To reduce the amount of bound charge on the charge terminalT₁, the desired elevation to provide free charge on the charge terminalT₁ for launching a guided surface wave should be at least 4-5 times theeffective spherical diameter above the lossy conductive medium (e.g.,the Earth). The compensation terminal T₂ can be used to adjust the totaleffective height (h_(TE)) of the guided surface waveguide probe 200 d toexcite an electric field having a guided surface wave tilt at R_(x). Thecompensation terminal T₂ can be positioned below the charge terminal T₁at h_(d)=h_(T)−h_(p), where h_(T) is the total physical height of thecharge terminal T₁. With the position of the compensation terminal T₂fixed and the phase delay Φ_(U) applied to the upper charge terminal T₁,the phase delay Φ_(L) applied to the lower compensation terminal T₂ canbe determined using the relationships of Equation (86), such that:

$\begin{matrix}{{\Phi_{U}\left( h_{d} \right)} = {{- {\beta \left( {h_{d} + y} \right)}} - {j\mspace{14mu} {{\ln\left( \frac{{R_{x} \times W} - {\left( {h_{p} + z} \right)e^{j{({{\beta \; h_{p}} + {\beta \; z} + \Phi_{L}})}}}}{\left( {h_{p} + y} \right)} \right)}.}}}} & (89)\end{matrix}$

In alternative embodiments, the compensation terminal T₂ can bepositioned at a height h_(d) where Im{Φ_(L)}=0. This is graphicallyillustrated in FIG. 15A, which shows plots 172 and 175 of the imaginaryand real parts of Φ_(U), respectively. The compensation terminal T₂ ispositioned at a height h_(d) where Im{Φ_(U)}=0, as graphicallyillustrated in plot 172. At this fixed height, the coil phase Φ_(U) canbe determined from Re{Φ_(U)}, as graphically illustrated in plot 175.

With the AC source 212 coupled to the coil 215 (e.g., at the 50Ω pointto maximize coupling), the position of tap 233 may be adjusted forparallel resonance of the compensation terminal T₂ with at least aportion of the coil at the frequency of operation. FIG. 15B shows aschematic diagram of the general electrical hookup of FIG. 14 in whichV₁ is the voltage applied to the lower portion of the coil 215 from theAC source 212 through tap 227, V₂ is the voltage at tap 224 that issupplied to the upper charge terminal T₁, and V₃ is the voltage appliedto the lower compensation terminal T₂ through tap 233. The resistancesR_(p) and R_(d) represent the ground return resistances of the chargeterminal T₁ and compensation terminal T₂, respectively. The charge andcompensation terminals T₁ and T₂ may be configured as spheres,cylinders, toroids, rings, hoods, or any other combination of capacitivestructures. The size of the charge and compensation terminals T₁ and T₂can be chosen to provide a sufficiently large surface for the charges Q₁and Q₂ imposed on the terminals. In general, it is desirable to make thecharge terminal T₁ as large as practical. The size of the chargeterminal T₁ should be large enough to avoid ionization of thesurrounding air, which can result in electrical discharge or sparkingaround the charge terminal. The self-capacitance C_(p) and C_(d) of thecharge and compensation terminals T₁ and T₂ respectively, can bedetermined using, for example, equation (24).

As can be seen in FIG. 15B, a resonant circuit is formed by at least aportion of the inductance of the coil 215, the self-capacitance C_(d) ofthe compensation terminal T₂, and the ground return resistance R_(d)associated with the compensation terminal T₂. The parallel resonance canbe established by adjusting the voltage V₃ applied to the compensationterminal T₂ (e.g., by adjusting a tap 233 position on the coil 215) orby adjusting the height and/or size of the compensation terminal T₂ toadjust C_(d). The position of the coil tap 233 can be adjusted forparallel resonance, which will result in the ground current through theground stake 218 and through the ammeter 236 reaching a maximum point.After parallel resonance of the compensation terminal T₂ has beenestablished, the position of the tap 227 for the AC source 212 can beadjusted to the 50Ω point on the coil 215.

Voltage V₂ from the coil 215 can be applied to the charge terminal T₁,and the position of tap 224 can be adjusted such that the phase (Φ) ofthe total effective height (h_(TE)) approximately equals the angle ofthe guided surface wave tilt (W_(Rx)) at the Hankel crossover distance(R_(x)). The position of the coil tap 224 can be adjusted until thisoperating point is reached, which results in the ground current throughthe ammeter 236 increasing to a maximum. At this point, the resultantfields excited by the guided surface waveguide probe 200 d aresubstantially mode-matched to a guided surface waveguide mode on thesurface of the lossy conducting medium 203, resulting in the launchingof a guided surface wave along the surface of the lossy conductingmedium 203. This can be verified by measuring field strength along aradial extending from the guided surface waveguide probe 200.

Resonance of the circuit including the compensation terminal T₂ maychange with the attachment of the charge terminal T₁ and/or withadjustment of the voltage applied to the charge terminal T₁ through tap224. While adjusting the compensation terminal circuit for resonanceaids the subsequent adjustment of the charge terminal connection, it isnot necessary to establish the guided surface wave tilt (W_(Rx)) at theHankel crossover distance (R_(x)). The system may be further adjusted toimprove coupling by iteratively adjusting the position of the tap 227for the AC source 212 to be at the 50Ω point on the coil 215 andadjusting the position of tap 233 to maximize the ground current throughthe ammeter 236. Resonance of the circuit including the compensationterminal T₂ may drift as the positions of taps 227 and 233 are adjusted,or when other components are attached to the coil 215.

In other implementations, the voltage V₂ from the coil 215 can beapplied to the charge terminal T₁, and the position of tap 233 can beadjusted such that the phase (Φ) of the total effective height (h_(TE))approximately equals the angle (Ψ) of the guided surface wave tilt atR_(x). The position of the coil tap 224 can be adjusted until theoperating point is reached, resulting in the ground current through theammeter 236 substantially reaching a maximum. The resultant fields aresubstantially mode-matched to a guided surface waveguide mode on thesurface of the lossy conducting medium 203, and a guided surface wave islaunched along the surface of the lossy conducting medium 203. This canbe verified by measuring field strength along a radial extending fromthe guided surface waveguide probe 200. The system may be furtheradjusted to improve coupling by iteratively adjusting the position ofthe tap 227 for the AC source 212 to be at the 50Ω point on the coil 215and adjusting the position of tap 224 and/or 233 to maximize the groundcurrent through the ammeter 236.

Referring back to FIG. 12, operation of a guided surface waveguide probe200 may be controlled to adjust for variations in operational conditionsassociated with the guided surface waveguide probe 200. For example, aprobe control system 230 can be used to control the feed network 209and/or positioning of the charge terminal T₁ and/or compensationterminal T₂ to control the operation of the guided surface waveguideprobe 200. Operational conditions can include, but are not limited to,variations in the characteristics of the lossy conducting medium 203(e.g., conductivity σ and relative permittivity ε_(r)), variations infield strength and/or variations in loading of the guided surfacewaveguide probe 200. As can be seen from Equations (41)-(44), the indexof refraction (n), the complex Brewster angle (θ_(i,B) and ψ_(i,B)), thewave tilt (|W|e^(jΨ)) and the complex effective height(h_(eff)=h_(p)e^(jΦ)) can be affected by changes in soil conductivityand permittivity resulting from, e.g., weather conditions.

Equipment such as, e.g., conductivity measurement probes, permittivitysensors, ground parameter meters, field meters, current monitors and/orload receivers can be used to monitor for changes in the operationalconditions and provide information about current operational conditionsto the probe control system 230. The probe control system 230 can thenmake one or more adjustments to the guided surface waveguide probe 200to maintain specified operational conditions for the guided surfacewaveguide probe 200. For instance, as the moisture and temperature vary,the conductivity of the soil will also vary. Conductivity measurementprobes and/or permittivity sensors may be located at multiple locationsaround the guided surface waveguide probe 200. Generally, it would bedesirable to monitor the conductivity and/or permittivity at or aboutthe Hankel crossover distance R_(x) for the operational frequency.Conductivity measurement probes and/or permittivity sensors may belocated at multiple locations (e.g., in each quadrant) around the guidedsurface waveguide probe 200.

With reference then to FIG. 16, shown is an example of a guided surfacewaveguide probe 200 e that includes a charge terminal T₁ and a chargeterminal T₂ that are arranged along a vertical axis z. The guidedsurface waveguide probe 200 e is disposed above a lossy conductingmedium 203, which makes up Region 1. In addition, a second medium 206shares a boundary interface with the lossy conducting medium 203 andmakes up Region 2. The charge terminals T₁ and T₂ are positioned overthe lossy conducting medium 203. The charge terminal T₁ is positioned atheight H₁, and the charge terminal T₂ is positioned directly below T₁along the vertical axis z at height H₂, where H₂ is less than H₁. Theheight h of the transmission structure presented by the guided surfacewaveguide probe 200 e is h=H₁−H₂. The guided surface waveguide probe 200e includes a feed network 209 that couples an excitation source 212 tothe charge terminals T₁ and T₂.

The charge terminals T₁ and/or T₂ include a conductive mass that canhold an electrical charge, which may be sized to hold as much charge aspractically possible. The charge terminal T₁ has a self-capacitance C₁,and the charge terminal T₂ has a self-capacitance C₂, which can bedetermined using, for example, equation (24). By virtue of the placementof the charge terminal T₁ directly above the charge terminal T₂, amutual capacitance C_(M) is created between the charge terminals T₁ andT₂. Note that the charge terminals T₁ and T₂ need not be identical, buteach can have a separate size and shape, and can include differentconducting materials. Ultimately, the field strength of a guided surfacewave launched by a guided surface waveguide probe 200 e is directlyproportional to the quantity of charge on the terminal T₁. The charge Q₁is, in turn, proportional to the self-capacitance C₁ associated with thecharge terminal T₁ since Q₁=C₁V, where V is the voltage imposed on thecharge terminal T₁.

When properly adjusted to operate at a predefined operating frequency,the guided surface waveguide probe 200 e generates a guided surface wavealong the surface of the lossy conducting medium 203. The excitationsource 212 can generate electrical energy at the predefined frequencythat is applied to the guided surface waveguide probe 200 e to excitethe structure. When the electromagnetic fields generated by the guidedsurface waveguide probe 200 e are substantially mode-matched with thelossy conducting medium 203, the electromagnetic fields substantiallysynthesize a wave front incident at a complex Brewster angle thatresults in little or no reflection. Thus, the surface waveguide probe200 e does not produce a radiated wave, but launches a guided surfacetraveling wave along the surface of a lossy conducting medium 203. Theenergy from the excitation source 212 can be transmitted as Zennecksurface currents to one or more receivers that are located within aneffective transmission range of the guided surface waveguide probe 200e.

One can determine asymptotes of the radial Zenneck surface currentJ_(ρ)(ρ) on the surface of the lossy conducting medium 203 to be J₁(ρ)close-in and J₂(ρ) far-out, where

$\begin{matrix}{{{{Close}\text{-}{in}\mspace{14mu} \left( {\rho < {\lambda \text{/}8}} \right)\text{:}\mspace{14mu} {J_{\rho}(\rho)}\text{∼}J_{1}} = {\frac{I_{1} + I_{2}}{2\; \pi \; \rho} + \frac{{E_{\rho}^{QS}\left( Q_{1} \right)} + {E_{\rho}^{QS}\left( Q_{2} \right)}}{Z_{\rho}}}},{and}} & (90) \\{\mspace{79mu} {{{Far}\text{-}{out}\mspace{14mu} \left( {\rho {\lambda \text{/}8}} \right)\text{:}\mspace{14mu} {J_{\rho}(\rho)}\text{∼}J_{2}} = {\frac{j\; \gamma \; \omega \; Q_{1}}{4} \times \sqrt{\frac{2\; \gamma}{\pi}} \times {\frac{e^{{- {({\alpha + {j\; \beta}})}}\rho}}{\sqrt{\rho}}.}}}} & (91)\end{matrix}$

where I₁ is the conduction current feeding the charge Q₁ on the firstcharge terminal T₁, and I₂ is the conduction current feeding the chargeQ₂ on the second charge terminal T₂. The charge Q₁ on the upper chargeterminal T₁ is determined by Q₁=C₁V₁, where C₁ is the isolatedcapacitance of the charge terminal T₁. Note that there is a thirdcomponent to J₁ set forth above given by (E_(ρ) ^(Q) ¹ )/Z_(ρ), whichfollows from the Leontovich boundary condition and is the radial currentcontribution in the lossy conducting medium 203 pumped by thequasi-static field of the elevated oscillating charge on the firstcharge terminal Q₁. The quantity Z_(ρ)=jωμ_(o)/γ_(e) is the radialimpedance of the lossy conducting medium, whereγ_(e)=(jωμ₁σ₁−ω²μ₁∈₁)^(1/2).

The asymptotes representing the radial current close-in and far-out asset forth by equations (90) and (91) are complex quantities. Accordingto various embodiments, a physical surface current J(ρ), is synthesizedto match as close as possible the current asymptotes in magnitude andphase. That is to say close-in, |J(ρ)| is to be tangent to |J₁|, andfar-out |J(ρ)| is to be tangent to |J₂|. Also, according to the variousembodiments, the phase of J(ρ) should transition from the phase of J₁close-in to the phase of J₂ far-out.

In order to match the guided surface wave mode at the site oftransmission to launch a guided surface wave, the phase of the surfacecurrent |J₂|far-out should differ from the phase of the surface current|J₁| close-in by the propagation phase corresponding to e^(−jβ(ρ) ²^(−ρ) ¹ ⁾ plus a constant of approximately 45 degrees or 225 degrees.This is because there are two roots for √{square root over (γ)}, onenear π/4 and one near 5π/4. The properly adjusted synthetic radialsurface current is

$\begin{matrix}{{J_{\rho}\left( {\rho,\varphi,0} \right)} = {\frac{I_{o}\gamma}{4}\mspace{14mu} {{H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}.}}} & (92)\end{matrix}$

Note that this is consistent with equation (17). By Maxwell's equations,such a J(ρ) surface current automatically creates fields that conform to

$\begin{matrix}{{H_{\varphi} = {\frac{{- \gamma}\; I_{o}}{4}e^{{- u_{2}}z}\mspace{14mu} {H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}}},} & (93) \\{{E_{\rho} = {\frac{{- \gamma}\; I_{o}}{4}\left( \frac{u_{2}}{j\; \omega \; ɛ_{0}} \right)e^{{- u_{2}}z}\mspace{14mu} {H_{1}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}}},{and}} & (94) \\{E_{z} = {\frac{{- \gamma}\; I_{o}}{4}\left( \frac{- \gamma}{\omega \; ɛ_{0}} \right)e^{{- u_{2}}z}\mspace{14mu} {{H_{0}^{(2)}\left( {{- j}\; \gamma \; \rho} \right)}.}}} & (95)\end{matrix}$

Thus, the difference in phase between the surface current |J₂| far-outand the surface current |J₁| close-in for the guided surface wave modethat is to be matched is due to the characteristics of the Hankelfunctions in equations (93)-(95), which are consistent with equations(1)-(3). It is of significance to recognize that the fields expressed byequations (1)-(6) and (17) and equations (92)-(95) have the nature of atransmission line mode bound to a lossy interface, not radiation fieldsthat are associated with groundwave propagation.

In order to obtain the appropriate voltage magnitudes and phases for agiven design of a guided surface waveguide probe 200 e at a givenlocation, an iterative approach may be used. Specifically, analysis maybe performed of a given excitation and configuration of a guided surfacewaveguide probe 200 e taking into account the feed currents to theterminals T₁ and T₂, the charges on the charge terminals T₁ and T₂, andtheir images in the lossy conducting medium 203 in order to determinethe radial surface current density generated. This process may beperformed iteratively until an optimal configuration and excitation fora given guided surface waveguide probe 200 e is determined based ondesired parameters. To aid in determining whether a given guided surfacewaveguide probe 200 e is operating at an optimal level, a guided fieldstrength curve 103 (FIG. 1) may be generated using equations (1)-(12)based on values for the conductivity of Region 1 (σ₁) and thepermittivity of Region 1 (ε₁) at the location of the guided surfacewaveguide probe 200 e. Such a guided field strength curve 103 canprovide a benchmark for operation such that measured field strengths canbe compared with the magnitudes indicated by the guided field strengthcurve 103 to determine if optimal transmission has been achieved.

In order to arrive at an optimized condition, various parametersassociated with the guided surface waveguide probe 200 e may beadjusted. One parameter that may be varied to adjust the guided surfacewaveguide probe 200 e is the height of one or both of the chargeterminals T₁ and/or T₂ relative to the surface of the lossy conductingmedium 203. In addition, the distance or spacing between the chargeterminals T₁ and T₂ may also be adjusted. In doing so, one may minimizeor otherwise alter the mutual capacitance C_(M) or any boundcapacitances between the charge terminals T₁ and T₂ and the lossyconducting medium 203 as can be appreciated. The size of the respectivecharge terminals T₁ and/or T₂ can also be adjusted. By changing the sizeof the charge terminals T₁ and/or T₂, one will alter the respectiveself-capacitances C₁ and/or C₂, and the mutual capacitance C_(M) as canbe appreciated.

Still further, another parameter that can be adjusted is the feednetwork 209 associated with the guided surface waveguide probe 200 e.This may be accomplished by adjusting the size of the inductive and/orcapacitive reactances that make up the feed network 209. For example,where such inductive reactances comprise coils, the number of turns onsuch coils may be adjusted. Ultimately, the adjustments to the feednetwork 209 can be made to alter the electrical length of the feednetwork 209, thereby affecting the voltage magnitudes and phases on thecharge terminals T₁ and T₂.

Note that the iterations of transmission performed by making the variousadjustments may be implemented by using computer models or by adjustingphysical structures as can be appreciated. By making the aboveadjustments, one can create corresponding “close-in” surface current J₁and “far-out” surface current J₂ that approximate the same currents J(ρ)of the guided surface wave mode specified in Equations (90) and (91) setforth above. In doing so, the resulting electromagnetic fields would besubstantially or approximately mode-matched to a guided surface wavemode on the surface of the lossy conducting medium 203.

While not shown in the example of FIG. 16, operation of the guidedsurface waveguide probe 200 e may be controlled to adjust for variationsin operational conditions associated with the guided surface waveguideprobe 200. For example, a probe control system 230 shown in FIG. 12 canbe used to control the feed network 209 and/or positioning and/or sizeof the charge terminals T₁ and/or T₂ to control the operation of theguided surface waveguide probe 200 e. Operational conditions caninclude, but are not limited to, variations in the characteristics ofthe lossy conducting medium 203 (e.g., conductivity σ and relativepermittivity ε_(r)), variations in field strength and/or variations inloading of the guided surface waveguide probe 200 e.

Referring now to FIG. 17, shown is an example of the guided surfacewaveguide probe 200 e of FIG. 16, denoted herein as guided surfacewaveguide probe 200 f. The guided surface waveguide probe 200 f includesthe charge terminals T₁ and T₂ that are positioned along a vertical axisz that is substantially normal to the plane presented by the lossyconducting medium 203 (e.g., the Earth). The second medium 206 is abovethe lossy conducting medium 203. The charge terminal T₁ has aself-capacitance C₁, and the charge terminal T₂ has a self-capacitanceC₂. During operation, charges Q₁ and Q₂ are imposed on the chargeterminals T₁ and T₂, respectively, depending on the voltages applied tothe charge terminals T₁ and T₂ at any given instant. A mutualcapacitance C_(M) may exist between the charge terminals T₁ and T₂depending on the distance there between. In addition, bound capacitancesmay exist between the respective charge terminals T₁ and T₂ and thelossy conducting medium 203 depending on the heights of the respectivecharge terminals T₁ and T₂ with respect to the lossy conducting medium203.

The guided surface waveguide probe 200 f includes a feed network 209that comprises an inductive impedance comprising a coil L_(1a) having apair of leads that are coupled to respective ones of the chargeterminals T₁ and T₂. In one embodiment, the coil L_(1a) is specified tohave an electrical length that is one-half (½) of the wavelength at theoperating frequency of the guided surface waveguide probe 200 f.

While the electrical length of the coil L_(1a) is specified asapproximately one-half (½) the wavelength at the operating frequency, itis understood that the coil L_(1a) may be specified with an electricallength at other values. According to one embodiment, the fact that thecoil L_(1a) has an electrical length of approximately one-half thewavelength at the operating frequency provides for an advantage in thata maximum voltage differential is created on the charge terminals T₁ andT₂. Nonetheless, the length or diameter of the coil L_(1a) may beincreased or decreased when adjusting the guided surface waveguide probe200 f to obtain optimal excitation of a guided surface wave mode.Adjustment of the coil length may be provided by taps located at one orboth ends of the coil. In other embodiments, it may be the case that theinductive impedance is specified to have an electrical length that issignificantly less than or greater than ½ the wavelength at theoperating frequency of the guided surface waveguide probe 200 f.

The excitation source 212 can be coupled to the feed network 209 by wayof magnetic coupling. Specifically, the excitation source 212 is coupledto a coil LP that is inductively coupled to the coil L_(1a). This may bedone by link coupling, a tapped coil, a variable reactance, or othercoupling approach as can be appreciated. To this end, the coil LP actsas a primary, and the coil L_(1a) acts as a secondary as can beappreciated.

In order to adjust the guided surface waveguide probe 200 f for thetransmission of a desired guided surface wave, the heights of therespective charge terminals T₁ and T₂ may be altered with respect to thelossy conducting medium 203 and with respect to each other. Also, thesizes of the charge terminals T₁ and T₂ may be altered. In addition, thesize of the coil L_(1a) may be altered by adding or eliminating turns orby changing some other dimension of the coil L_(1a). The coil L_(1a) canalso include one or more taps for adjusting the electrical length asshown in FIG. 17. The position of a tap connected to either chargeterminal T₁ or T₂ can also be adjusted.

Referring next to FIGS. 18A, 18B, 18C and 19, shown are examples ofgeneralized receive circuits for using the surface-guided waves inwireless power delivery systems. FIGS. 18A and 18B-18C include a linearprobe 303 and a tuned resonator 306, respectively. FIG. 19 is a magneticcoil 309 according to various embodiments of the present disclosure.According to various embodiments, each one of the linear probe 303, thetuned resonator 306, and the magnetic coil 309 may be employed toreceive power transmitted in the form of a guided surface wave on thesurface of a lossy conducting medium 203 according to variousembodiments. As mentioned above, in one embodiment the lossy conductingmedium 203 comprises a terrestrial medium (or Earth).

With specific reference to FIG. 18A, the open-circuit terminal voltageat the output terminals 312 of the linear probe 303 depends upon theeffective height of the linear probe 303. To this end, the terminalpoint voltage may be calculated as

V _(T)=∫₀ ^(h) ^(e) E _(inc) ·dl,  (96)

where E_(inc) is the strength of the incident electric field induced onthe linear probe 303 in Volts per meter, dl is an element of integrationalong the direction of the linear probe 303, and h_(e) is the effectiveheight of the linear probe 303. An electrical load 315 is coupled to theoutput terminals 312 through an impedance matching network 318.

When the linear probe 303 is subjected to a guided surface wave asdescribed above, a voltage is developed across the output terminals 312that may be applied to the electrical load 315 through a conjugateimpedance matching network 318 as the case may be. In order tofacilitate the flow of power to the electrical load 315, the electricalload 315 should be substantially impedance matched to the linear probe303 as will be described below.

Referring to FIG. 18B, a ground current excited coil 306 a possessing aphase shift equal to the wave tilt of the guided surface wave includes acharge terminal T_(R) that is elevated (or suspended) above the lossyconducting medium 203. The charge terminal T_(R) has a self-capacitanceC_(R). In addition, there may also be a bound capacitance (not shown)between the charge terminal T_(R) and the lossy conducting medium 203depending on the height of the charge terminal T_(R) above the lossyconducting medium 203. The bound capacitance should preferably beminimized as much as is practicable, although this may not be entirelynecessary in every instance.

The tuned resonator 306 a also includes a receiver network comprising acoil L_(R) having a phase shift Φ. One end of the coil L_(R) is coupledto the charge terminal T_(R), and the other end of the coil L_(R) iscoupled to the lossy conducting medium 203. The receiver network caninclude a vertical supply line conductor that couples the coil L_(R) tothe charge terminal T_(R). To this end, the coil L_(R) (which may alsobe referred to as tuned resonator L_(R)-C_(R)) comprises aseries-adjusted resonator as the charge terminal C_(R) and the coilL_(R) are situated in series. The phase delay of the coil L_(R) can beadjusted by changing the size and/or height of the charge terminalT_(R), and/or adjusting the size of the coil L_(R) so that the phase Φof the structure is made substantially equal to the angle of the wavetilt Ψ. The phase delay of the vertical supply line can also be adjustedby, e.g., changing length of the conductor.

For example, the reactance presented by the self-capacitance C_(R) iscalculated as 1/jωC_(R). Note that the total capacitance of thestructure 306 a may also include capacitance between the charge terminalT_(R) and the lossy conducting medium 203, where the total capacitanceof the structure 306 a may be calculated from both the self-capacitanceC_(R) and any bound capacitance as can be appreciated. According to oneembodiment, the charge terminal T_(R) may be raised to a height so as tosubstantially reduce or eliminate any bound capacitance. The existenceof a bound capacitance may be determined from capacitance measurementsbetween the charge terminal T_(R) and the lossy conducting medium 203 aspreviously discussed.

The inductive reactance presented by a discrete-element coil L_(R) maybe calculated as jωL, where L is the lumped-element inductance of thecoil L_(R). If the coil L_(R) is a distributed element, its equivalentterminal-point inductive reactance may be determined by conventionalapproaches. To tune the structure 306 a, one would make adjustments sothat the phase delay is equal to the wave tilt for the purpose ofmode-matching to the surface waveguide at the frequency of operation.Under this condition, the receiving structure may be considered to be“mode-matched” with the surface waveguide. A transformer link around thestructure and/or an impedance matching network 324 may be insertedbetween the probe and the electrical load 327 in order to couple powerto the load. Inserting the impedance matching network 324 between theprobe terminals 321 and the electrical load 327 can effect aconjugate-match condition for maximum power transfer to the electricalload 327.

When placed in the presence of surface currents at the operatingfrequencies power will be delivered from the surface guided wave to theelectrical load 327. To this end, an electrical load 327 may be coupledto the structure 306 a by way of magnetic coupling, capacitive coupling,or conductive (direct tap) coupling. The elements of the couplingnetwork may be lumped components or distributed elements as can beappreciated.

In the embodiment shown in FIG. 18B, magnetic coupling is employed wherea coil L_(S) is positioned as a secondary relative to the coil L_(R)that acts as a transformer primary. The coil L_(S) may be link-coupledto the coil L_(R) by geometrically winding it around the same corestructure and adjusting the coupled magnetic flux as can be appreciated.In addition, while the receiving structure 306 a comprises aseries-tuned resonator, a parallel-tuned resonator or even adistributed-element resonator of the appropriate phase delay may also beused.

While a receiving structure immersed in an electromagnetic field maycouple energy from the field, it can be appreciated thatpolarization-matched structures work best by maximizing the coupling,and conventional rules for probe-coupling to waveguide modes should beobserved. For example, a TE₂₀ (transverse electric mode) waveguide probemay be optimal for extracting energy from a conventional waveguideexcited in the TE₂₀ mode. Similarly, in these cases, a mode-matched andphase-matched receiving structure can be optimized for coupling powerfrom a surface-guided wave. The guided surface wave excited by a guidedsurface waveguide probe 200 on the surface of the lossy conductingmedium 203 can be considered a waveguide mode of an open waveguide.Excluding waveguide losses, the source energy can be completelyrecovered. Useful receiving structures may be E-field coupled, H-fieldcoupled, or surface-current excited.

The receiving structure can be adjusted to increase or maximize couplingwith the guided surface wave based upon the local characteristics of thelossy conducting medium 203 in the vicinity of the receiving structure.To accomplish this, the phase delay (Φ) of the receiving structure canbe adjusted to match the angle (Ψ) of the wave tilt of the surfacetraveling wave at the receiving structure. If configured appropriately,the receiving structure may then be tuned for resonance with respect tothe perfectly conducting image ground plane at complex depth z=−d/2.

For example, consider a receiving structure comprising the tunedresonator 306 a of FIG. 18B, including a coil L_(R) and a verticalsupply line connected between the coil L_(R) and a charge terminalT_(R). With the charge terminal T_(R) positioned at a defined heightabove the lossy conducting medium 203, the total phase shift Φ of thecoil L_(R) and vertical supply line can be matched with the angle (Ψ) ofthe wave tilt at the location of the tuned resonator 306 a. FromEquation (22), it can be seen that the wave tilt asymptotically passesto

$\begin{matrix}{{W = {{{W}e^{j\; \Psi}} = {\frac{E_{\rho}}{E_{z}}\underset{\rho\rightarrow\infty}{\rightarrow}\frac{1}{\sqrt{ɛ_{r} - {j\frac{\sigma_{1}}{\omega \; ɛ_{o}}}}}}}},} & (97)\end{matrix}$

where ε_(r) comprises the relative permittivity and σ₁ is theconductivity of the lossy conducting medium 203 at the location of thereceiving structure, ε_(o) is the permittivity of free space, and ω=2πf,where f is the frequency of excitation. Thus, the wave tilt angle (Ψ)can be determined from Equation (97).

The total phase shift (Φ=θ_(c)+θ_(y)) of the tuned resonator 306 aincludes both the phase delay (θ_(c)) through the coil L_(R) and thephase delay of the vertical supply line (θ_(y)). The spatial phase delayalong the conductor length l_(w) of the vertical supply line can begiven by θ_(y)=β_(w)l_(w), where β_(w) is the propagation phase constantfor the vertical supply line conductor. The phase delay due to the coil(or helical delay line) is θ_(c)=β_(p)l_(C), with a physical length ofl_(C) and a propagation factor of

$\begin{matrix}{{\beta_{p} = {\frac{2\; \pi}{\lambda_{p}} = \frac{2\; \pi}{V_{f}\lambda_{0}}}},} & (98)\end{matrix}$

where V is the velocity factor on the structure, λ₀ is the wavelength atthe supplied frequency, and λ_(p) is the propagation wavelengthresulting from the velocity factor V_(f). One or both of the phasedelays (θ_(c)+θ_(y)) can be adjusted to match the phase shift Φ to theangle (Ψ) of the wave tilt. For example, a tap position may be adjustedon the coil L_(R) of FIG. 18B to adjust the coil phase delay (θ_(c)) tomatch the total phase shift to the wave tilt angle (Φ=Ψ). For example, aportion of the coil can be bypassed by the tap connection as illustratedin FIG. 18B. The vertical supply line conductor can also be connected tothe coil L_(R) via a tap, whose position on the coil may be adjusted tomatch the total phase shift to the angle of the wave tilt.

Once the phase delay (Φ) of the tuned resonator 306 a has been adjusted,the impedance of the charge terminal T_(R) can then be adjusted to tuneto resonance with respect to the perfectly conducting image ground planeat complex depth z=−d/2. This can be accomplished by adjusting thecapacitance of the charge terminal T₁ without changing the travelingwave phase delays of the coil L_(R) and vertical supply line. Theadjustments are similar to those described with respect to FIGS. 9A and9B.

The impedance seen “looking down” into the lossy conducting medium 203to the complex image plane is given by:

Z _(in) =R _(in) +jX _(in) =Z _(o) tan h(jβ _(o)(d/2)),  (99)

where β_(o)=ω√{square root over (μ_(o)∈_(o))}. For vertically polarizedsources over the Earth, the depth of the complex image plane can begiven by:

d/2≈1/√{square root over (jωμ ₁σ₁−ω²μ₁ε₁)},  (100)

where μ₁ is the permeability of the lossy conducting medium 203 andε₁=ε_(r)ε_(o).

At the base of the tuned resonator 306 a, the impedance seen “lookingup” into the receiving structure is Z_(↑)=Z_(base) as illustrated inFIG. 9A. With a terminal impedance of:

$\begin{matrix}{{Z_{R} = \frac{1}{j\; \omega \; C_{R}}},} & (101)\end{matrix}$

where C_(R) is the self-capacitance of the charge terminal T_(R), theimpedance seen “looking up” into the vertical supply line conductor ofthe tuned resonator 306 a is given by:

$\begin{matrix}{{Z_{2} = {{Z_{W}\frac{Z_{R} + {Z_{w}{\tanh \left( {j\; \beta_{w}h_{w}} \right)}}}{Z_{w} + {Z_{R}{\tanh \left( {j\; \beta_{w}h_{w}} \right)}}}} = {Z_{W}\frac{Z_{R} + {Z_{w}{\tanh \left( {j\; \theta_{y}} \right)}}}{Z_{w} + {Z_{R}{\tanh \left( {j\; \theta_{y}} \right)}}}}}},} & (102)\end{matrix}$

and the impedance seen “looking up” into the coil L_(R) of the tunedresonator 306 a is given by:

$\begin{matrix}{Z_{base} = {{R_{base} + {jX}_{base}} = {{Z_{R}\frac{Z_{2} + {Z_{R}{\tanh \left( {j\; \beta_{p}H} \right)}}}{Z_{R} + {Z_{2}\tanh \; \left( {j\; \beta_{p}H} \right)}}} = {Z_{c}{\frac{Z_{2} + {Z_{R}{\tanh \left( {j\; \theta_{c}} \right)}}}{Z_{R} + {Z_{2}{\tanh \left( {j\; \theta_{c}} \right)}}}.}}}}} & (103)\end{matrix}$

By matching the reactive component (X_(in)) seen “looking down” into thelossy conducting medium 203 with the reactive component (X_(base)) seen“looking up” into the tuned resonator 306 a, the coupling into theguided surface waveguide mode may be maximized.

Referring next to FIG. 18C, shown is an example of a tuned resonator 306b that does not include a charge terminal T_(R) at the top of thereceiving structure. In this embodiment, the tuned resonator 306 b doesnot include a vertical supply line coupled between the coil L_(R) andthe charge terminal T_(R). Thus, the total phase shift (Φ) of the tunedresonator 306 b includes only the phase delay (θ_(c)) through the coilL_(R). As with the tuned resonator 306 a of FIG. 18B, the coil phasedelay θ_(c) can be adjusted to match the angle (Ψ) of the wave tiltdetermined from Equation (97), which results in Φ=Ψ. While powerextraction is possible with the receiving structure coupled into thesurface waveguide mode, it is difficult to adjust the receivingstructure to maximize coupling with the guided surface wave without thevariable reactive load provided by the charge terminal T_(R).

Referring to FIG. 18D, shown is a flow chart 180 illustrating an exampleof adjusting a receiving structure to substantially mode-match to aguided surface waveguide mode on the surface of the lossy conductingmedium 203. Beginning with 181, if the receiving structure includes acharge terminal T_(R) (e.g., of the tuned resonator 306 a of FIG. 18B),then the charge terminal T_(R) is positioned at a defined height above alossy conducting medium 203 at 184. As the surface guided wave has beenestablished by a guided surface waveguide probe 200, the physical height(h_(p)) of the charge terminal T_(R) may be below that of the effectiveheight. The physical height may be selected to reduce or minimize thebound charge on the charge terminal T_(R) (e.g., four times thespherical diameter of the charge terminal). If the receiving structuredoes not include a charge terminal T_(R) (e.g., of the tuned resonator306 b of FIG. 18C), then the flow proceeds to 187.

At 187, the electrical phase delay Φ of the receiving structure ismatched to the complex wave tilt angle Ψ defined by the localcharacteristics of the lossy conducting medium 203. The phase delay(θ_(c)) of the helical coil and/or the phase delay (θ_(y)) of thevertical supply line can be adjusted to make Φ equal to the angle (Ψ) ofthe wave tilt (W). The angle (Ψ) of the wave tilt can be determined fromEquation (86). The electrical phase Φ can then be matched to the angleof the wave tilt. For example, the electrical phase delay Φ=θ_(c)+θ_(y)can be adjusted by varying the geometrical parameters of the coil L_(R)and/or the length (or height) of the vertical supply line conductor.

Next at 190, the load impedance of the charge terminal T_(R) can betuned to resonate the equivalent image plane model of the tunedresonator 306 a. The depth (d/2) of the conducting image ground plane139 (FIG. 9A) below the receiving structure can be determined usingEquation (100) and the values of the lossy conducting medium 203 (e.g.,the Earth) at the receiving structure, which can be locally measured.Using that complex depth, the phase shift (θ_(d)) between the imageground plane 139 and the physical boundary 136 (FIG. 9A) of the lossyconducting medium 203 can be determined using θ_(d)=β₀ d/2. Theimpedance (Z_(in)) as seen “looking down” into the lossy conductingmedium 203 can then be determined using Equation (99). This resonancerelationship can be considered to maximize coupling with the guidedsurface waves.

Based upon the adjusted parameters of the coil L_(R) and the length ofthe vertical supply line conductor, the velocity factor, phase delay,and impedance of the coil L_(R) and vertical supply line can bedetermined. In addition, the self-capacitance (C_(R)) of the chargeterminal T_(R) can be determined using, e.g., Equation (24). Thepropagation factor (β_(p)) of the coil L_(R) can be determined usingEquation (98), and the propagation phase constant (β_(w)) for thevertical supply line can be determined using Equation (49). Using theself-capacitance and the determined values of the coil L_(R) andvertical supply line, the impedance (Z_(base)) of the tuned resonator306 a as seen “looking up” into the coil L_(R) can be determined usingEquations (101), (102), and (103).

The equivalent image plane model of FIG. 9A also applies to the tunedresonator 306 a of FIG. 18B. The tuned resonator 306 a can be tuned toresonance with respect to the complex image plane by adjusting the loadimpedance Z_(R) of the charge terminal T_(R) such that the reactancecomponent X_(base) of Z_(base) cancels out the reactance component ofX_(in) of Z_(in), or X_(base)+X_(in)=0. Thus, the impedance at thephysical boundary 136 (FIG. 9A) “looking up” into the coil of the tunedresonator 306 a is the conjugate of the impedance at the physicalboundary 136 “looking down” into the lossy conducting medium 203. Theload impedance Z_(R) can be adjusted by varying the capacitance (C_(R))of the charge terminal T_(R) without changing the electrical phase delayΦ=θ_(c)+θ_(y) seen by the charge terminal T_(R). An iterative approachmay be taken to tune the load impedance Z_(R) for resonance of theequivalent image plane model with respect to the conducting image groundplane 139. In this way, the coupling of the electric field to a guidedsurface waveguide mode along the surface of the lossy conducting medium203 (e.g., Earth) can be improved and/or maximized.

Referring to FIG. 19, the magnetic coil 309 comprises a receive circuitthat is coupled through an impedance matching network 333 to anelectrical load 336. In order to facilitate reception and/or extractionof electrical power from a guided surface wave, the magnetic coil 309may be positioned so that the magnetic flux of the guided surface wave,H_(φ), passes through the magnetic coil 309, thereby inducing a currentin the magnetic coil 309 and producing a terminal point voltage at itsoutput terminals 330. The magnetic flux of the guided surface wavecoupled to a single turn coil is expressed by

=∫∫_(A) _(CS) μ_(r)μ_(o) {right arrow over (H)}·{circumflex over(n)}dA  (104)

where

is the coupled magnetic flux, μ_(r) is the effective relativepermeability of the core of the magnetic coil 309, μ_(o) is thepermeability of free space, {right arrow over (H)} is the incidentmagnetic field strength vector, {circumflex over (n)} is a unit vectornormal to the cross-sectional area of the turns, and A_(CS) is the areaenclosed by each loop. For an N-turn magnetic coil 309 oriented formaximum coupling to an incident magnetic field that is uniform over thecross-sectional area of the magnetic coil 309, the open-circuit inducedvoltage appearing at the output terminals 330 of the magnetic coil 309is

$\begin{matrix}{{V = {{{- N}\frac{d\; \mathcal{F}}{dt}} \approx {{- j}\; \omega \; \mu_{r}\mu_{0}{NHA}_{CS}}}},} & (105)\end{matrix}$

where the variables are defined above. The magnetic coil 309 may betuned to the guided surface wave frequency either as a distributedresonator or with an external capacitor across its output terminals 330,as the case may be, and then impedance-matched to an external electricalload 336 through a conjugate impedance matching network 333.

Assuming that the resulting circuit presented by the magnetic coil 309and the electrical load 336 are properly adjusted and conjugateimpedance matched, via impedance matching network 333, then the currentinduced in the magnetic coil 309 may be employed to optimally power theelectrical load 336. The receive circuit presented by the magnetic coil309 provides an advantage in that it does not have to be physicallyconnected to the ground.

With reference to FIGS. 18A, 18B, 18C and 19, the receive circuitspresented by the linear probe 303, the mode-matched structure 306, andthe magnetic coil 309 each facilitate receiving electrical powertransmitted from any one of the embodiments of guided surface waveguideprobes 200 described above. To this end, the energy received may be usedto supply power to an electrical load 315/327/336 via a conjugatematching network as can be appreciated. This contrasts with the signalsthat may be received in a receiver that were transmitted in the form ofa radiated electromagnetic field. Such signals have very low availablepower, and receivers of such signals do not load the transmitters.

It is also characteristic of the present guided surface waves generatedusing the guided surface waveguide probes 200 described above that thereceive circuits presented by the linear probe 303, the mode-matchedstructure 306, and the magnetic coil 309 will load the excitation source212 (e.g., FIGS. 3, 12 and 16) that is applied to the guided surfacewaveguide probe 200, thereby generating the guided surface wave to whichsuch receive circuits are subjected. This reflects the fact that theguided surface wave generated by a given guided surface waveguide probe200 described above comprises a transmission line mode. By way ofcontrast, a power source that drives a radiating antenna that generatesa radiated electromagnetic wave is not loaded by the receivers,regardless of the number of receivers employed.

Thus, together one or more guided surface waveguide probes 200 and oneor more receive circuits in the form of the linear probe 303, the tunedmode-matched structure 306, and/or the magnetic coil 309 can make up awireless distribution system. Given that the distance of transmission ofa guided surface wave using a guided surface waveguide probe 200 as setforth above depends upon the frequency, it is possible that wirelesspower distribution can be achieved across wide areas and even globally.

The conventional wireless-power transmission/distribution systemsextensively investigated today include “energy harvesting” fromradiation fields and also sensor coupling to inductive or reactivenear-fields. In contrast, the present wireless-power system does notwaste power in the form of radiation which, if not intercepted, is lostforever. Nor is the presently disclosed wireless-power system limited toextremely short ranges as with conventional mutual-reactance couplednear-field systems. The wireless-power system disclosed hereinprobe-couples to the novel surface-guided transmission line mode, whichis equivalent to delivering power to a load by a waveguide or a loaddirectly wired to the distant power generator. Not counting the powerrequired to maintain transmission field strength plus that dissipated inthe surface waveguide, which at extremely low frequencies isinsignificant relative to the transmission losses in conventionalhigh-tension power lines at 60 Hz, all of the generator power goes onlyto the desired electrical load. When the electrical load demand isterminated, the source power generation is relatively idle.

Referring next to FIGS. 20A-E, shown are examples of various schematicsymbols that are used with reference to the discussion that follows.With specific reference to FIG. 20A, shown is a symbol that representsany one of the guided surface waveguide probes 200 a, 200 b, 200 c, 200e, 200 d, or 200 f, or any variations thereof. In the following drawingsand discussion, a depiction of this symbol will be referred to as aguided surface waveguide probe P. For the sake of simplicity in thefollowing discussion, any reference to the guided surface waveguideprobe P is a reference to any one of the guided surface waveguide probes200 a, 200 b, 200 c, 200 e, 200 d, or 200 f, or variations thereof.

Similarly, with reference to FIG. 20B, shown is a symbol that representsa guided surface wave receive structure that may comprise any one of thelinear probe 303 (FIG. 18A), the tuned resonator 306 (FIGS. 18B-18C), orthe magnetic coil 309 (FIG. 19). In the following drawings anddiscussion, a depiction of this symbol will be referred to as a guidedsurface wave receive structure R. For the sake of simplicity in thefollowing discussion, any reference to the guided surface wave receivestructure R is a reference to any one of the linear probe 303, the tunedresonator 306, or the magnetic coil 309; or variations thereof.

Further, with reference to FIG. 20C, shown is a symbol that specificallyrepresents the linear probe 303 (FIG. 18A). In the following drawingsand discussion, a depiction of this symbol will be referred to as aguided surface wave receive structure R_(P). For the sake of simplicityin the following discussion, any reference to the guided surface wavereceive structure R_(P) is a reference to the linear probe 303 orvariations thereof.

Further, with reference to FIG. 20D, shown is a symbol that specificallyrepresents the tuned resonator 306 (FIGS. 18B-18C). In the followingdrawings and discussion, a depiction of this symbol will be referred toas a guided surface wave receive structure R_(R). For the sake ofsimplicity in the following discussion, any reference to the guidedsurface wave receive structure R_(R) is a reference to the tunedresonator 306 or variations thereof.

Further, with reference to FIG. 20E, shown is a symbol that specificallyrepresents the magnetic coil 309 (FIG. 19). In the following drawingsand discussion, a depiction of this symbol will be referred to as aguided surface wave receive structure R_(M). For the sake of simplicityin the following discussion, any reference to the guided surface wavereceive structure R_(M) is a reference to the magnetic coil 309 orvariations thereof.

With reference to FIG. 21, shown is a global transmission environment400. The global transmission environment 400 includes a computingenvironment 403, one or more emergency identification devices 406, aguided surface waveguide probe P, a disaster warning device 412, anetwork 415, and a guided surface wave 418.

According to one embodiment, the computing environment 403 includes oneor more data store(s) 421, an emergency application 424, and a powermodulator 425. The data store 421 includes one or more data elements,such as remote locations 427, which may include, among other data,coordinates 430, language 433, and name 436. The disaster warning device412 includes guided surface wave receive structure R, power circuitry443, data processing circuitry 446, and warning output circuitry 449.The power modulator 425 includes an excitation source and ultimatelyexperiences the disaster warning device 412 as a load. The powermodulator 425 may be a power supply with a circuit configured tomodulate a frequency carrier waveform containing the data onto the powertransmission to send a message as will be described.

The one or more emergency identification devices 406 and the emergencyapplication 424 communicate via the network 415. The one or moreemergency identification devices 406 identify that a global disasterevent is presently occurring or has occurred, or may predict a disasterevent will likely occur. The one or more emergency identificationdevices 406 identify the global disaster event by being notified from atrusted third-party source, through measuring one or more environmentalsensing or prediction devices or through data entry from a user. The oneor more emergency identification devices 406 indicate the occurrence ofthe global disaster event by transmitting a message to the emergencyapplication 424 via the network 415. The message includes details aboutone or more global disaster events including location informationdescribing where the one or more global disaster events taking place andthe location that may be affected by the global disaster event.

For example, the message may include details of an earthquake presentlyoccurring or that recently occurred in the Pacific Ocean nearCalifornia, USA with coordinates defining the area of the earthquake.The message may include coordinates defining an area of concern alongthe Pacific Rim based on a potential risk of tsunami. Such a message mayinclude details of when the tsunami may occur. The global disasterevents may include, among other disasters, one or more earthquakes,tsunamis, hurricanes, tornadoes, solar flares, volcanic eruptions,lightning storms, landslides, cyclones, typhoons, floods, limniceruptions, blizzards, droughts, hailstorms, heat waves, wildfires,impact from space such as an asteroid, a bomb threat, a terroristthreat, and health disasters such as a disease outbreak and/or epidemic.

In one embodiment, the emergency application 424 authenticates andvalidates the global disaster event. For example, the emergencyapplication 424 may require that an indication of a global disasterevent be transmitted using a secure protocol, such as SSL, andauthenticated via a third party, such as a trusted certificateauthority. In some embodiments, the network 415 is a secure and/orprivate network. The emergency application 424 may require that theindication be confirmed by more than one of the one or more emergencyidentification devices 406, as discussed herein.

In another embodiment, the emergency application 424 directs the powermodulator 425 to perform one or more actions, such as by communicatingdata over a guided surface wave 418 via the guided surface waveguideprobe P. The one or more actions may include transmitting and receivingdata via guided surface wave 418 on one or more specified, negotiated,or predetermined frequencies or transmitting power via guided surfacewave 418 on one or more specified or predetermined frequencies. The oneor more actions may include ceasing to transmit power via guided surfacewaves 418 on one or more specified, negotiated, or predeterminedfrequencies. The data may be transmitted or received by the powermodulator 425 via the guided surface waveguide probe P by modulating orembedding a data signal into the power signal of the guided surface wave418.

The disaster warning device 412 receives power and data from one or moreguided surface waveguide probes P by way of one or more guided surfacewaves 418 received by the guided surface wave receive structure R. Insome embodiments, the power and the data are combined into a singleguided surface wave 418 using an appropriate modulation scheme. Theguided surface wave receive structure R is electrically coupled to thepower circuitry 443 and the data processing circuitry 446. The powercircuitry 443 may supply power to the disaster warning device 412. Insome embodiments, the power circuitry 443 obtains power from a combinedpower signal and data signal transmitted, for example, by the guidedsurface waveguide probe P.

According to one embodiment, the data processing circuitry 446 receivesdata from a guided surface wave 418. For example, the data processingcircuitry 446 extracts a modulated data signal embedded in a guidedsurface wave 418 that includes a combined power signal and data signal.The data processing circuitry 446 is electrically coupled to the warningoutput circuitry 449. The warning output circuitry 449 includes, amongother output devices, one or more auditory devices, such as a speaker,and/or one or more visual displays, such as an LCD display or a beaconlight.

According to one embodiment, the data processing circuitry 446demodulates and validates information about one or more global disasterevents, and outputs a warning to users proximate to the disaster warningdevice 412 using the warning output circuitry 449 via the warning outputcircuitry 449. In some embodiments, the disaster warning device 412 isconfigured to warn users within a proximate distance of several hundredfeet from the device. In other embodiments, the disaster warning device412 may include warning output circuitry 449 capable of outputting awarning message within a proximate distance of miles, or tens of milesbased on the geographical requirements associated with a given area andthe likelihood and severity of a potential global disaster event in aparticular area.

With reference to FIGS. 22-24, shown are embodiments of the disasterwarning device 412, denoted herein as disaster warning devices 412A and412B, according to various embodiments. The circuits shown are merelyexamples of various embodiments of the present disclosure. Differentdesigns may also be used to accomplish the same functionality.

In FIG. 22, according to one embodiment, the disaster warning device412A includes guided surface wave receive structure R, power circuitry443, data processing circuitry 446A, and warning output circuitry 449.The disaster warning device 412A also includes a positive outputterminal 543, a negative output terminal 546, and an electrical coupling549. In some embodiments, the data processing circuitry 446A includes anoperational amplifier 529, a signal processing circuitry 533, and acomputing device 536. The computing device 536 is executes a disasterapplication 537. The operational amplifier 529 includes an outputconnection 550, a positive power connection 552, a negative powerconnection 555, a positive input connection 558, and a negative inputconnection 561. The output 550 of the operational amplifier 529 isconnected to the signal processing circuitry 533 and connected to ground523 via a resistor.

According to one embodiment, the power circuitry 443 includes animpedance matching network 509, a power processing circuitry 512, and apower storage circuitry 513. The impedance matching network 509 may bethe same as any of impedance matching networks 318, 324, and 333 andminimizes reflections and provides for maximum power transfer as can beappreciated. The power storage circuitry 513 includes one or morecapacitors and/or one or more batteries. In some embodiments, the powerstorage circuitry 513 stores power from the power processing circuitry512, such as a DC power signal for use by the data processing circuitry446A and the warning output circuitry 449. The power processingcircuitry 512 includes a transformer, a voltage rectifier, and/or avoltage regulator. The power circuitry 443 converts an AC signal fromthe guided surface wave receive structure R into a DC circuit. To thisend, the power circuitry 443 may include a power rectifier.

According to one embodiment, the power storage circuitry 513 includes abattery to supply power to the disaster warning device 412 for aspecific period of time. The specified period of time may exceed apredetermined threshold time period that a guided surface waveguideprobe P is configured to transmit a guided surface wave on a frequencymatching the guide surface wave receive structure R. For example, thepower storage circuitry 513 may be configured to supply power to thedisaster warning device 412, such as the disaster warning device 412A or412B, for at least 24 hours, and the guided surface waveguide probe Psupplies enough power to charge the power storage circuitry 513 every 6hours, thereby ensuring the disaster warning device 412 does not run outof power during normal operation. The guided surface waveguide probe Pmay continuously power the disaster warning device 412 with the powerstorage circuitry 513 storing power as a backup.

In some embodiments, the power circuitry 443 includes circuitryconfigured to decide whether or not to draw power based on one or morefactors. The power circuitry 443 draws power only during a preconfiguredinterval or window of time. The interval or window of time is configuredvia a message transmitted to data processing circuitry 446 or dataprocessing circuitry 515 (FIG. 23). For example, the disasterapplication 537 receives a message indicating a new interval or windowof time to draw power and cause the power circuitry to be configured forthe new interval or window of time. The disaster application 537 mayconfigure the new interval or window of time by setting a memoryreadable by the power circuitry 443. In one embodiment, the new intervalor window of time is a time period to return the power storage circuitry513 to a fully charged state. The charging of the power storagecircuitry 513 to the fully charged state may be done periodically at atime offset with the charging of batteries of other devices. The numberof disaster warning device 412A charging at one time may be based inpart on a capacity of power available. The power circuitry 443 mayinclude a computing device communicably coupled to the computing device536, such a computing device in the data processing circuitry 515.

According to one embodiment, the operational amplifier 529 provides aninfinite resistance on the positive input connection 558. Theoperational amplifier 529 functions as a buffer for output of the guidedsurface wave receive structure R, wherein the output 550 may mirror thepositive input connection 558 without disturbing the output of theguided surface wave receive structure R when processed by the signalprocessing circuitry 533. The positive input connection 558 is connectedto an output of the guided surface wave receive structure R, as shown inFIGS. 22 and 24, or is connected to an output of the power circuitry443, such as output 543, as shown in FIG. 23.

The signal processing circuitry 533 is electrically coupled to theoutput 550 of the operational amplifier 529. According to oneembodiment, the signal processing circuitry 533 is configured todemodulate data from a guided surface wave received via the guidedsurface wave receive structure R, such as the guided surface wave 418(FIG. 21). To do so, the signal processing circuitry 533 includes anAD/DC converter to convert the demodulated data signal from an analogsignal to a digital data signal from which information can be obtained.The signal processing circuitry 533 processes the digital data, forexample, converting the digital data into data packets, and transmitsthe data to a disaster application 537. The disaster application 537processes the incoming data and extract information describing one ormore global disaster events. The disaster application 537 extracts orobtains geographic information describing one or more locations and/orareas that are affected by the one or more global disaster events andcompare the geographic information to localized geographic informationdescribing an area associated with the disaster warning device 412A.

For example, the disaster warning device 412A may be installed on abeach in Japan along the Eastern coast and have localized geographicinformation including the longitude and latitude of the installationlocation of the disaster warning device 412A. In this example, thedisaster warning device 412A identifies a message indicating that aglobal disaster event in Japan is intended for the disaster warningdevice 412A by determining that the longitude and latitude of thedisaster warning device 412A are within a geographic area described inthe geographic information from the message indicating the globaldisaster event. Further, the disaster warning device 412A identifies amessage identifying a global disaster event in France as not intendedfor the disaster warning device 412A by determining the longitude andlatitude of the disaster warning device 412A are outside of thegeographic area. In one embodiment, the disaster warning device 412Aincludes a GPS device configured to identify the installation locationof the disaster warning device 412A.

According to one embodiment, the computing device 536 outputsinformation about a global disaster event using the warning outputcircuitry 449. The warning output circuitry includes an output device550. The output device 550 may include one or more of a speaker, an airhorn, a light, an LCD display, a billboard, and/or other device. Theoutput device 550 may initiate a safety procedure and/or initiate asafety device. A safety device may include raising a driving barrier,activating a switch, raising or lowering a bridge, launching a drone,opening a dam, transmitting an electromagnetic pulse, and/or launchingan offensive weapon. For example, in response to receiving a warningthat a tsunami will approach a beach in 30 minutes, the disaster warningdevice 412 may initiate converting all lanes for roads nearby to adirection leading away from the beach via an electric lane controlsystem to minimize traffic, flash warning lights on the beach shore,output an audio message to evacuate the beach, and launch a drone toprovide video support for recovery specialists, such as police,firefighters, medical workers, and the military.

According to one embodiment, the disaster warning device 412 includes anamplification of the input signal to improve the incoming signal priorto processing. In this embodiment, the operational amplifier 529includes a first resistor connected between the negative inputconnection 561 and the output 550, and a second resistor connectedbetween the output 550 and the ground 523. In some embodiments, thesignal at the output 550 may be equal to the positive input connection558 multiplied by the result of the first resistor divided by the secondresistor. In some embodiments, a value of the first resistor may beprogrammatically adjusted to tune the signal output from the operationalamplifier 529 based on the input range of the signal processingcircuitry 533. The disaster application 537 may be configured to adjustthe value of the first resistor based on the output of the signalprocessing circuitry 533 and data describing the signal level receivedby the signal processing circuitry 533 from the operational amplifier529.

The disaster warning device 412A senses and/or receives a guided surfacewave on guided surface wave receive structure R. The guided surface wavereceive structure R outputs a signal as a result. Once the signal isgenerated, the power circuitry 443 generates power from the signal. Thepower is used to power the disaster warning device 412A in real time orby storing the generated power signal in the power storage circuitry513. In some embodiments, the disaster warning device 412A buffers theincoming signal using a buffer circuit, such as the operationalamplifier 529. The guided surface wave signal may be buffered before orafter the power circuitry 443. In some embodiments, the buffer circuitmay amplify the signal to improve the strength of the signal.

As shown in FIG. 23, the disaster warning device 412B may include morethan one guided surface wave receive structure R, shown as guidedsurface wave receive structures R1 and R2. The guided surface wavereceive structures R1 and R2 are each a guided surface wave receivestructures R configured to operate on a frequency different from eachother. The guide surface wave receive structure R1 may provide a signalfor the power circuitry 443. The data processing circuitry 446 mayreceive a signal on a frequency differing from the frequency used topower the disaster warning device 412B. The guide surface wave receivestructure R2 may receive a guided surface wave and output the signal onthe connection 573 to the data processing circuitry 446A.

According to one embodiment, the power processing circuitry 512 includesa data processing circuitry 515 communicably connected to the guidedsurface wave receive structure R1 configured to receive data via aguided surface wave using the frequency of the guided surface wavereceive structure R1. The data processing circuitry 515 may be a dataprocessing circuitry 446, such as data processing circuitry 446A or446B. Data communications via the guided surface wave receive structureR1 may be limited to power characteristics of the power circuitry 512.The data processing circuitry 515 may transmit power usagecharacteristics to the computing environment 403 via a wireless or wirednetwork. For example, the second data processing circuitry 446 mayreport a usage of 1 KWH over a period of one month.

The data processing circuitry 515 may set a new interval or window oftime during which the power circuitry 443B draws power from the powersignal received on the guided surface wave receive structure R1. Thedata processing circuitry 515 may receive a message disabling thedisaster warning device 412B. In some embodiments, the power circuitry512 may include circuitry to draw enough power while disabled to operatethe data processing circuitry 515. In this embodiment, while disabled,the disaster warning device 412B becomes fully enabled in response to amessage from the computing environment 403.

According to one embodiment, the disaster warning device 412B receives aguided surface wave 418 on a first frequency predetermined fortransmission of power via the guided surface wave receive structure R1.The disaster warning device 412B generates a power signal usingelectrical energy from the guided surface wave receive structure R1. Thedisaster warning device 412B is experienced as a load at an excitationsource coupled to the guided surface waveguide probe P (FIG. 21). Insome embodiments, the power modulator 425 (FIG. 21) includes and/or isthe excitation source.

According to one embodiment, the power generation from the powercircuitry 443B is experienced as a load at power modulator 425. Thepower processing circuitry 512 stores the generated power in the powerstorage circuitry 513 for future use. The disaster warning device 412Bmay receive a guided surface wave 418 on a second frequencypredetermined for the transmission of data via the guided surface wavereceive structure R2.

The power storage circuitry 513 provides power to various elements ofthe disaster warning device 412B, such as the data processing circuitry446A and the warning output 449. In some embodiments, the dataprocessing circuitry 515 communicates via a wireless network withemergency application 424. The data processing circuitry 515 receivescommands from emergency application 424 via one or more of the guidedsurface wave receive structure R1 or the guided surface wave receivestructure R2 to regulate power storage and power use.

The emergency application 424 directs the data processing circuitry 515to store and provide differing amounts of power, which can be based onnumerous factors including current power demands on guided surfacewaveguide probe P, time of day, the occurrence of a disaster event, aquantity of power required, and other factors. In one embodiment, theemergency application 424 directs the data processing circuitry 515, bymodulating a data signal via power modulator 425, to refrain fromdrawing power for a period of time in response to a disaster eventoccurring in an area determined not to include the data processingcircuitry 515. The disaster warning device 412B uses power from thepower storage circuitry 513 to power data processing circuitry 446 andthe warning output 449 in order to process the second guided surfacewave.

Referring next to FIG. 24, shown is a flowchart that provides oneexample of the operation of a portion of the emergency application 424according to various embodiments. It is understood that the flowchart ofFIG. 24 provides merely an example of the many different types offunctional arrangements that may be employed to implement the operationof the portion of the emergency application 424 as described herein. Asan alternative, the flowchart of FIG. 24 may be viewed as depicting anexample of elements of a method implemented in the computing environment403 (FIG. 21) according to one or more embodiments.

The following flowchart illustrates how the emergency application 424receives a message indicating an emergency event from an authenticatedsource. The emergency application 424 determines one or more locationsaffected by the emergency event and generates emergency messages withcoordinates of locations affected. The emergency application 424transmits the generated emergency messages and receives anacknowledgement from the disaster warning devices 412.

At box 576, the emergency application 424 initiates the process. At box579, the emergency application 424 includes receiving a messageindicating that an emergency event has occurred. For example, theemergency application 424 may receive a message from the emergencyidentification device 406 via the network 415 indicating that anemergency has occurred. The message may include a type of emergencyevent, a location of the emergency event, one or more areas of concernfor future emergency events, areas of risk for future emergency events,a risk of occurrence for one or more future emergency events, a severitylevel of the emergency event, or a predicted severity level for a futureemergency event. The message may also include security information, suchas a digital signature and/or encryption data.

At box 582, the emergency application 424 includes authenticating thesource of the message indicating the emergency event. For example, theemergency identification device 406 may have a public security key, suchas an RSS Key, signed by a trusted third party. The emergencyapplication 424 may authenticate the emergency identification device 406by ensuring that the message can be decrypted by a public key that isdigitally signed by the trusted third party. In some embodiments, theoccurrence of the emergency event may be authenticated by a process suchas the one described in FIG. 26.

At box 585, the emergency application 424 includes determining thelocations affected by an emergency event. For example, the authenticatedmessage from the emergency identification device 406 may includeinformation describing affected areas or locations. The emergencyapplication 424 may extract the information describing the affectedareas or locations. In some embodiments, the emergency application 424searches the data store 421 for the one or more remote locations 427that have the coordinates 430 within the information describing theaffected areas or locations. The remote locations 427 may comprise thelocations of the one or more disaster warning devices 412 including acoordinates 430, a language 433, and a name 436 of associated with thecorresponding disaster warning device 412.

At box 588, the emergency application 424 includes generating one ormore emergency messages including coordinates for affected locations.For example, the emergency application 424 may generate a message foreach of the determined locations in remote locations 427 including datadescribing the emergency event. The emergency application 424 may groupthe remote locations 427, such as the disaster warning devices 412,together based on one or more characteristic of the location. Theemergency application 424 creates groups of remote locations in whichgroup members are chosen for a group based on the group members beingwithin a certain distance from one another or within a preconfigureddistance from the emergency event. The emergency application 424 alsocreates groups of remote locations based on all members of the groupsharing a language 433. The emergency application 424 may generatemultiple messages individually corresponding to each remote location 427or each group of remote locations 427. In some embodiments, the messageswill include a description of the emergency event in the languagematching the one or more remote locations 427 to which the message isintended to be sent.

At box 591, the emergency application 424 includes transmitting anemergency message via a guided surface waveguide probe P. For example,the emergency application 424 transmits an emergency message via aguided surface waveguide probe P (FIG. 21) by causing the guided surfacewaveguide probe P to modulate a data signal onto a guided surface waveat a predetermined frequency, such as via power modulator 425 (FIG. 21).The emergency application 424 directs the guided surface waveguide probeP via a device driver or may communicate the message to a computingdevice electrically coupled to the guided surface waveguide probe P. Thefrequency of transmission is specified so that the corresponding guidedsurface wave transmitted will reach the respective disaster warningdevice(s) 412. In one embodiment, the frequency is specified to be lowenough to facilitate global propagation of the guided surface wave.

At box 594, the emergency application 424 may include determining if oneor more of the remote devices received the guided surface wave. Forexample, the emergency application 424 may receive a message via network415 indicating that a remote device has acknowledged receiving a messagevia a guided surface wave. If the message is not acknowledged, theemergency application 424 may retransmit the message at box 591. In someembodiments, the message is retransmitted a set number of times beforeproceeding to box 597 regardless. If a remote device fails toacknowledge a message after a predefined threshold of failed retransmitsor if a predefined threshold of messages fail to be acknowledged, theemergency application 424 may designate the remote device as inactive.The emergency application 424 may also transmit a notification when aremote device is designated as inactive. The emergency application 424may also present a user interface detailing the current status and anycurrent error conditions for all remote devices. If the message isacknowledged, the emergency application 424 ends at box 597.

Referring next to FIG. 25, shown is a flowchart that provides oneexample of the operation of a portion of the emergency application 424according to various embodiments. It is understood that the flowchart ofFIG. 25 provides merely an example of the many different types offunctional arrangements that may be employed to implement the operationof the portion of the emergency application 424 as described herein. Asan alternative, the flowchart of FIG. 25 may be viewed as depicting anexample of elements of a method implemented in the computing environment403 (FIG. 21) according to one or more embodiments.

The following flowchart describes how the emergency application 424determines locations affected by an emergency event and looks up one ormore codes describing areas within the affected area. To provideextensive coverage while minimizing data, the emergency application 424categorizes and/or scores each of the one or more codes. Then, theemergency application 424 minimizes a quantity of the one or more codesbased on the category and/or score.

At box 600, the emergency application 424 initiates the process. At box602, the emergency application 424 involves determining locationsaffected by an emergency event. For example, the emergency application424 may perform as described in box 585 in FIG. 24.

At box 604, the emergency application 424 involves determining one ormore codes describing areas within the affected area. According to oneembodiment, the emergency application 424 compares the affected area toa list of locations or areas in remote locations 427 (FIG. 21). Theremote locations 427 may be individually mapped to individual disasterwarning devices 412. Some remote locations 427 may cover groups ofdisaster warning devices 412. A single disaster warning device 412 maybe responsive to numerous codes. For example, a disaster warning device412 on a coast of a country may correspond to a code for the deviceitself, a code for the coast line, a code for areas surrounding anocean, and a code for the country itself. In addition, the disasterwarning device 412 may correspond to codes for smaller areas within thecoast line. The areas affected by respective potential disasters maydictate which disaster warning devices 412 are responsive to whichcodes. For example, one code may correspond to a potential tsunami onthe West coast of the United States. As such, all disaster warningdevices 412 along the West coast of the United States may be responsiveto this code.

At box 606, the emergency application 424 involves categorizing and/orscoring each of the one or more codes. For example, the emergencyapplication 424 may categorize, and/or read the categorization from datastore 421 (FIG. 21), each of the determined one or more codes describingareas within the affected area. The categories and/or scores may bebased on how much of the affected area is within a coverage area for acorresponding code. For example, a code corresponding to a coastal rimmay be less appropriate for use to describe a fire along a small area ofthe coast, even if a small subset of the coastal rim is affected.

The categories and/or scores may include and/or be based in part on asubset of the coverage area for a code is affected and/or the entirecoverage area is affected. The categories and/or scores may also includeand/or be based in part on the subset of the affected area that thecoverage area for a code covers. For example, hundreds of disasterwarning devices 412 may be included in an affected area, and a code thatcorresponds to 1, 5, or 10 of the disaster warning devices 412 is scoredlower than a code that corresponds to 100 or more disaster warningdevices 412. The scores may also be based in part on an interrelation oftwo or more codes relative to the affected area. For example, twodifferent codes may each receive a higher score if, when combined, thecoverage area encompasses a substantial part of the affected area or allof the affected area.

At box 608, the emergency application 424 involves minimizing thequantity of one or more codes based on the categories and/or scores. Theemergency application 424 identifies a most efficient subset of the oneor more codes that covers all disaster warning devices within theaffected area. For example, the emergency application 424 may combinecodes from within multiple categories to ensure all disaster warningdevices 412 within the affected area are notified. The emergencyapplication 424 may identify a code representing a group of devices butomitting one or more disaster warning devices 412 within the affectedarea and combine the code for the group with individual codes for theone or more disaster warning devices 412 omitted from the group.

The emergency application 424 may choose codes based on minimizing thenumber of disaster warning devices 412 notified that are outside of theaffected area. A number of disaster warning devices 412 may be includedin a group where a substantial portion of the group is located in theaffected area, but the number of disaster warning devices 412 are notwithin the affected area. The emergency application 424 may avoidselecting the code for that group if the number of disaster warningdevices 412 is large and/or select the code for that group to transmitif the number of disaster warning devices 412 is small relative to thenumber of disaster warning devices 412 in the affected area. Theemergency application 424 ends at box 610.

Referring next to FIG. 26, shown is a flowchart that provides oneexample of the operation of a portion of the disaster application 537(FIG. 22) according to various embodiments. It is understood that theflowchart of FIG. 26 provides merely an example of the many differenttypes of functional arrangements that may be employed to implement theoperation of the portion of the disaster application 537 as describedherein. As an alternative, the flowchart of FIG. 26 may be viewed asdepicting an example of elements of a method implemented in disasterapplication 537 according to one or more embodiments.

The following flowchart describes how the disaster application 537,which is executed within a disaster warning device 412, demodulates adata signal from a guided surface wave. Using the data signal, thedisaster application 537 determines an area affected by an emergencyevent. Next, the disaster application 537 checks to see if the disasterwarning device 412 is located within the affected area. If the disasterwarning device 412 is located within the affected area, a warning isoutput.

At box 612, the disaster application 537 initiates the process. At box615, the disaster application 537 includes demodulating a data signalfrom the guided surface wave signal. For example, the disasterapplication 537 may demodulate a data signal from the guided surfacewave signal using the signal processing circuitry 533. The signalprocessing circuitry 533 may provide the data to the disasterapplication 537 as a digital signal. The signal processing circuitry 533may convert the data signal into digital data and transmit it via theelectrical connection 540. In some embodiments, the electricalconnection 540 connects to a GPIO pin of the computing device 536, andthe computing device 536 converts the low level input into packetizedbytes or other format for processing by the disaster application 537.The disaster application 537 may acknowledge a message by transmittingan acknowledgement that the message was successfully received via awireless connection.

At box 618, the disaster application 537 includes determining anaffected area from the demodulated data signal. For example, thedisaster application 537 may process packetized bytes from the signalprocessing circuitry to extract data packets describing an emergencyevent. The data packets may include details of one or more areas orlocations affected by the emergency event. At box 621, the disasterapplication 537 includes determining if a location of the disasterwarning device 412 itself is in the affected area. The computing device536 may include a GPS device to determine its present coordinates andthe disaster application 537 may compare the present coordinates tocoordinates of an area of coordinates within the message.

Alternatively, the message may include addresses/unique identifierscorresponding to affected disaster warning devices 412, as described inFIGS. 25 and 27. If the device location is within the affected area(e.g., Yes is answered), the disaster application 537 may proceed to box624. If the device is not within the affected area (e.g. No isanswered), the disaster application 537 may proceed to box 627.

At box 624, the disaster application 537 includes outputting a warningto the sounded area of the device is located in the affected area. Forexample, the disaster application 537 may output a warning from thewarning output 449 using connection 549, such as from the one or moreoutput devices 550. For example, if a disaster warning device 412 werelocated on a beach where a Tsunami is about to hit, a message from theemergency application 424 may direct the disaster warning device 412 tooutput a warning message to people swimming at the beach instructingthem to leave the proximity of the beach because a Tsunami is coming. Insome embodiments, the output device 550 includes multiple auditoryand/or visual indicators covering a large area. For example, the beachmay have a speaker output, warning lights, and other items every 100yards along a coast wherein one or more disaster warning devices canoutput a warning on all of the speakers, lights, and other items. Uponsuccessfully outputting the warning message, the process proceeds to box627. It is possible that local devices may communicate locally (via awired or wireless channel) with a respective disaster warning device 412to further distribute the creation of a warning around the area.

Referring next to FIG. 27, shown is a flowchart that provides oneexample of the operation of a portion of the disaster application 537according to various embodiments. It is understood that the flowchart ofFIG. 27 provides merely an example of the many different types offunctional arrangements that may be employed to implement the operationof the portion of the disaster application 537 as described herein. Asan alternative, the flowchart of FIG. 27 may be viewed as depicting anexample of elements of a method implemented in the disaster warningdevice 412 (FIG. 21) according to one or more embodiments.

The following flowchart describes how the disaster application 537,which is executed within a disaster warning device 412, receives one ormore codes describing affected areas of an emergency event. The disasterapplication 537 also determines one or more codes assigned to thedisaster warning device 412. Next, the disaster application 537 comparesthe received one or more codes to the assigned one or more codes todetermine if the disaster warning device 412 is located within theaffected area. If the disaster warning device 412 is located within theaffected area, a warning is output.

Beginning at box 630, the disaster application 537 initiates theprocess. In some embodiments, the process occurs within box 618 of FIG.26. At box 633, the disaster application 537 involves receiving one ormore codes describing the affected areas of the emergency event. Thedisaster application 537 may identify one or more codes in a messagedescribing a disaster event.

At box 636, the disaster application 537 determines one or more codesapplicable to the disaster warning device 412. To this end, the disasterwarning device 412 may be responsive to one or more codes. One of theone or more codes may correspond only to the disaster warning 412.Another of the one or more codes may correspond to an area of which thedisaster warning device 412 is located. In some embodiments, thedisaster application 537 loads from memory one or more codes to whichthe disaster warning device 412 is responsive. In another embodiment,the disaster application 537 receives the one or more codes from aremote server via a wireless connection.

At box 639, the disaster application 537 determines if the receivedcodes match the codes assigned to the disaster warning device 412. Ifthe codes match, the disaster application 537 proceeds to box 642 (e.g.Yes). Otherwise (e.g. No), the disaster application 537 proceeds to box645. At box 642, the disaster application 537 outputs the warning, asdescribed in box 624 of FIG. 26.

Referring next to FIG. 28, shown is a flowchart that provides oneexample of the operation of a portion of the emergency application 424according to various embodiments. It is understood that the flowchart ofFIG. 28 provides merely an example of the many different types offunctional arrangements that may be employed to implement the operationof the portion of the emergency application 424 as described herein. Asan alternative, the flowchart of FIG. 28 may be viewed as depicting anexample of elements of a method implemented in the computing environment403 (FIG. 21) according to one or more embodiments.

FIG. 28 shows an example embodiment of the emergency application 424determining if a remoted emergency event actually occurred. Uponreceiving an initial indication that an emergency event has occurredfrom an external server 406, the emergency application 424 calculates anaccuracy threshold and queries other external servers 406 to confirm thereported emergency event is not fraudulent. If a calculated likelihoodof the event occurring exceeds the accuracy threshold, the emergencyapplication 424 transmits the emergency event message.

At box 654, the emergency application 424 initiates the process. At box657, the emergency application 424 involves receiving a first messageindicating an emergency event has occurred from a first source. Forexample, the emergency application 424 may receive a message from afirst emergency identification device 406 (FIG. 21) indicating that anemergency event has occurred. The message may include one or moredetails describing a disaster event including details of the severity ofthe event, areas affected by the event, the type of event, informationabout any current and potential damage in the affected area, and/ortiming of the event, such as already occurred, occurring now, will occurwith a percentage likelihood within a window, and/or will occur with apercentage likelihood at a calculated time.

At box 660, the emergency application 424 calculates an accuracythreshold. For example, the emergency application 424 may calculate anaccuracy threshold based in part on a number of factors includinginformation received in the message. The accuracy threshold may be lowerwhen, among other reasons, the risk of the disaster event is high, theseverity of the disaster event is high, and/or the time until occurrenceof the disaster event is short or has already occurred. The accuracythreshold may be higher when, among other reasons, the risk of thedisaster event is low, the severity of the disaster event is low, and/orthe time until occurrence of the disaster event is long.

At box 663, the emergency application 424 involves querying one or moresecond sources on the likelihood of occurrence of the reported emergencyevent. For example, the emergency application 424 may query one or moresecond emergency identification devices 406 requesting a likelihood thatthe disaster event occurred. The emergency application 424 may query theone or more second emergency identification devices 406 requesting toreport and/or analyze various data, such as seismographic data, weatherdata, news data, astrology data, ocean current data, and/or volcanodata, among other sources. As a non-limiting example, in response to areport of an incoming Tsunami from a first emergency identificationdevice 406, the emergency application 424 may query one or more secondemergency identification devices 406 requesting seismographic data tosee if an earthquake occurred along a coastline or in the ocean. Thequery may be a search request, such as a GOOGLE search request.

At box 666, the emergency application 424 involves receiving one or moresecond messages indicating a likelihood that the emergency eventoccurred from the one or more second sources. For example, the emergencyapplication 424 may receive one or more second messages from one or moreemergency identification devices 406. The second messages may be sent inresponse to a request, such as the request in box 663 and/or in responseto the emergency identification device 406 determining a disaster eventhas occurred. The emergency application 424 may receive the secondmessages resembling the first message from the second sources whenmultiple emergency identification devices 406 have detected the samedisaster event.

At box 669, the emergency application 424 involves determining alikelihood of occurrence that the emergency event occurred or willoccur. The determined likelihood may be based in part on the likelihoodreported by the one or more second sources and/or on informationobtained from the one or more second sources. The determined likelihoodmay be based on the credibility of the first emergency identificationdevice 406 and/or the credibility of the one or more second sources. Thecredibility may be calculated based on a history of interacting with asource and/or a predefined credibility.

At box 672, the emergency application 424 involves determining whetherthe likelihood of occurrence meets or exceeds the accuracy threshold.The emergency application 424 may determine that the event has occurredin response to determining that a determined likelihood meets theaccuracy threshold. In some embodiments, the emergency application 424may determine that an event has occurred in response to receiving morethan one messages indicating the same emergency event from the multiplesources. If the threshold is met, the emergency application 424 proceedsto box 675. Otherwise, the emergency application 424 proceeds to box678.

At box 675, the emergency application 424 involves initiatingtransmission of an emergency event message. For example, the emergencyapplication 424 may initiate a process including boxes 585, 588, 591,594, and 597 of FIG. 24. The emergency application 424 may initiate theprocess in response to determining that the likelihood of the eventoccurring meets or exceeds the accuracy threshold. Next, emergencyapplication 424 proceeds to box 678.

With reference to FIG. 29, shown is a schematic block diagram of acomputing device 681 in the global transmission environment 400 (FIG.21) according to an embodiment of the present disclosure. The globaltransmission environment 400 includes one or more computing environments403 and one or more computing devices 536. Each computing device 681includes at least one processor circuit, for example, having a processor684 and a memory 696, both of which are coupled to a local interface693. To this end, each computing device 681 may comprise, for example,at least one server computer or like device. The local interface 693 maycomprise, for example, a data bus with an accompanying address/controlbus or other bus structure as can be appreciated.

Stored in the memory 696 are both data and several components that areexecutable by the processor 684. In particular, stored in the memory 696and executable by the processor 684 are emergency application 424 (FIG.21) and disaster application 537 (FIG. 22), and potentially otherapplications. Also stored in the memory 696 may be a data store 421(FIG. 21) and other data. In addition, an operating system may be storedin the memory 696 and executable by the processor 684.

It is understood that there may be other applications that are stored inthe memory 696 and are executable by the processor 684 as can beappreciated. Where any component discussed herein is implemented in theform of software, any one of a number of programming languages may beemployed such as, for example, C, C++, C#, Objective C, Java®,JavaScript®, Perl, PHP, Visual Basic®, Python®, Ruby, Flash®, or otherprogramming languages.

A number of software components are stored in the memory 696 and areexecutable by the processor 684. In this respect, the term “executable”means a program file that is in a form that can ultimately be run by theprocessor 684. Examples of executable programs may be, for example, acompiled program that can be translated into machine code in a formatthat can be loaded into a random access portion of the memory 696 andrun by the processor 684, source code that may be expressed in properformat such as object code that is capable of being loaded into a randomaccess portion of the memory 696 and executed by the processor 684, orsource code that may be interpreted by another executable program togenerate instructions in a random access portion of the memory 696 to beexecuted by the processor 684, etc. An executable program may be storedin any portion or component of the memory 696 including, for example,random access memory (RAM), read-only memory (ROM), hard drive,solid-state drive, USB flash drive, memory card, optical disc such ascompact disc (CD) or digital versatile disc (DVD), floppy disk, magnetictape, or other memory components.

The memory 696 is defined herein as including both volatile andnonvolatile memory and data storage components. Volatile components arethose that do not retain data values upon loss of power. Nonvolatilecomponents are those that retain data upon a loss of power. Thus, thememory 696 may comprise, for example, random access memory (RAM),read-only memory (ROM), hard disk drives, solid-state drives, USB flashdrives, memory cards accessed via a memory card reader, floppy disksaccessed via an associated floppy disk drive, optical discs accessed viaan optical disc drive, magnetic tapes accessed via an appropriate tapedrive, and/or other memory components, or a combination of any two ormore of these memory components. In addition, the RAM may comprise, forexample, static random access memory (SRAM), dynamic random accessmemory (DRAM), or magnetic random access memory (MRAM) and other suchdevices. The ROM may comprise, for example, a programmable read-onlymemory (PROM), an erasable programmable read-only memory (EPROM), anelectrically erasable programmable read-only memory (EEPROM), or otherlike memory device.

Also, the processor 684 may represent multiple processors 684 and/ormultiple processor cores and the memory 696 may represent multiplememories 696 that operate in parallel processing circuits, respectively.In such a case, the local interface 693 may be an appropriate networkthat facilitates communication between any two of the multipleprocessors 684, between any processor 684 and any of the memories 696,or between any two of the memories 696, etc. The local interface 693 maycomprise additional systems designed to coordinate this communication,including, for example, performing load balancing. The processor 684 maybe of electrical or of some other available construction.

Although emergency application 424 and disaster application 537, andother various systems described herein may be embodied in software orcode executed by general purpose hardware as discussed above, as analternative the same may also be embodied in dedicated hardware or acombination of software/general purpose hardware and dedicated hardware.If embodied in dedicated hardware, each can be implemented as a circuitor state machine that employs any one of or a combination of a number oftechnologies. These technologies may include, but are not limited to,discrete logic circuits having logic gates for implementing variouslogic functions upon an application of one or more data signals,application specific integrated circuits (ASICs) having appropriatelogic gates, field-programmable gate arrays (FPGAs), or othercomponents, etc. Such technologies are generally well known by thoseskilled in the art and, consequently, are not described in detailherein.

The flowcharts of FIGS. 24, 25, and 28 show the functionality andoperation of an implementation of portions of the emergency application424 and the flowchart of FIGS. 26 and 27 show the functionality andoperation for an implementation of the disaster application 537. Ifembodied in software, each block may represent a module, segment, orportion of code that comprises program instructions to implement thespecified logical function(s). The program instructions may be embodiedin the form of source code that comprises human-readable statementswritten in a programming language or machine code that comprisesnumerical instructions recognizable by a suitable execution system suchas a processor 684 in a computer system or other system. The machinecode may be converted from the source code, etc. If embodied inhardware, each block may represent a circuit or a number ofinterconnected circuits to implement the specified logical function(s).

Although the flowcharts of FIGS. 24-28 show a specific order ofexecution, it is understood that the order of execution may differ fromthat which is depicted. For example, the order of execution of two ormore blocks may be scrambled relative to the order shown. Also, two ormore blocks shown in succession in FIGS. 24-28 may be executedconcurrently or with partial concurrence. Further, in some embodiments,one or more of the blocks shown in FIGS. 24-28 may be skipped oromitted. In addition, any number of counters, state variables, warningsemaphores, or messages might be added to the logical flow describedherein, for purposes of enhanced utility, accounting, performancemeasurement, or providing troubleshooting aids, etc. It is understoodthat all such variations are within the scope of the present disclosure.

Also, any logic or application described herein, including emergencyapplication 424 and disaster application 537, that comprises software orcode can be embodied in any non-transitory computer-readable medium foruse by or in connection with an instruction execution system such as,for example, a processor 684 in a computer system or other system. Inthis sense, the logic may comprise, for example, statements includinginstructions and declarations that can be fetched from thecomputer-readable medium and executed by the instruction executionsystem. In the context of the present disclosure, a “computer-readablemedium” can be any medium that can contain, store, or maintain the logicor application described herein for use by or in connection with theinstruction execution system.

The computer-readable medium can comprise any one of many physical mediasuch as, for example, magnetic, optical, or semiconductor media. Morespecific examples of a suitable computer-readable medium would include,but are not limited to, magnetic tapes, magnetic floppy diskettes,magnetic hard drives, memory cards, solid-state drives, USB flashdrives, or optical discs. Also, the computer-readable medium may be arandom access memory (RAM) including, for example, static random accessmemory (SRAM) and dynamic random access memory (DRAM), or magneticrandom access memory (MRAM). In addition, the computer-readable mediummay be a read-only memory (ROM), a programmable read-only memory (PROM),an erasable programmable read-only memory (EPROM), an electricallyerasable programmable read-only memory (EEPROM), or other type of memorydevice.

Further, any logic or application described herein, including emergencyapplication 424 and disaster application 537, may be implemented andstructured in a variety of ways. For example, one or more applicationsdescribed may be implemented as modules or components of a singleapplication. Further, one or more applications described herein may beexecuted in shared or separate computing devices or a combinationthereof. For example, a plurality of the applications described hereinmay execute in the same computing device 681, or in multiple computingdevices in the same global transmission environment 400. Additionally,it is understood that terms such as “application,” “service,” “system,”“engine,” “module,” and so on may be interchangeable and are notintended to be limiting.

Disjunctive language such as the phrase “at least one of X, Y, or Z,”unless specifically stated otherwise, is otherwise understood with thecontext as used in general to present that an item, term, etc., may beeither X, Y, or Z, or any combination thereof (e.g., X, Y, and/or Z).Thus, such disjunctive language is not generally intended to, and shouldnot, imply that certain embodiments require at least one of X, at leastone of Y, or at least one of Z to each be present.

It should be emphasized that the above-described embodiments of thepresent disclosure are merely possible examples of implementations setforth for a clear understanding of the principles of the disclosure.Many variations and modifications may be made to the above-describedembodiment(s) without departing substantially from the spirit andprinciples of the disclosure. All such modifications and variations areintended to be included herein within the scope of this disclosure andprotected by the following claims. In addition, all optional andpreferred features and modifications of the described embodiments anddependent claims are usable in all aspects of the disclosure taughtherein. Furthermore, the individual features of the dependent claims, aswell as all optional and preferred features and modifications of thedescribed embodiments are combinable and interchangeable with oneanother.

Therefore, the following is claimed:
 1. A system, comprising: a datastore; a power modulator coupled to a guided surface waveguide probe;and at least one computing device configured to at least: determine alocation affected by an emergency event; generate an emergency messagecorresponding to the emergency event; and transmit the emergency messagevia a guided surface wave launched by the guided surface waveguideprobe.
 2. The system of claim 1, wherein the power modulator isconfigured to embed a data signal into a power signal of the guidedsurface wave.
 3. The system of claim 1, wherein the power modulator isconfigured to excite a charge terminal of the guided surface waveguideprobe to generate an electromagnetic field that synthesizes a wave frontincident at a complex Brewster angle of incidence of a terrestrialmedium to generate the guided surface wave that travels across theterrestrial medium.
 4. The system of claim 3, further comprising a feednetwork electrically coupled between the charge terminal and the powermodulator, the feed network providing a phase delay that matches a wavetilt angle associated with the complex Brewster angle of incidence. 5.The system of claim 1, wherein the at least one computing device isfurther configured to at least receive a message indicating theemergency event occurred from a third party service.
 6. The system ofclaim 5, wherein the at least one computing device is further configuredto at least authenticate the third party service based at least in parton a digital signature.
 7. The system of claim 1, wherein the at leastone computing device is further configured to at least receive anacknowledgement from at least one disaster warning device.
 8. A system,comprising: a data store; and at least one computing device configuredto at least: receive a message indicating an emergency event occurredfrom a third party service; generate an emergency message correspondingto the emergency event; and transmit the emergency message via a guidedsurface wave launched by a guided surface waveguide probe.
 9. The systemof claim 8, wherein the at least one computing device is furtherconfigured to at least determine an affected area for the emergencyevent, wherein the emergency message comprises coordinates correspondingto the affected area.
 10. The system of claim 8, wherein the at leastone computing device is further configured to at least: determine anaffected area for the emergency event; and determine at least one codecorresponding to the affected area, wherein the emergency message isbased at least in part on the at least one code.
 11. The system of claim10, wherein the at least one code comprises a plurality of codes and theat least one computing device is further configured to at least:generate a respective score for each of the plurality of codes; andminimize a quantity of the plurality of codes based at least on part onthe respective score for each of the plurality of codes.
 12. The systemof claim 10, wherein the at least one code comprises a plurality ofcodes and the at least one computing device is further configured to atleast: categorize the plurality of codes; and minimize a quantity of theplurality of codes based at least on part on the categorization.
 13. Amethod, comprising: determining, via at least one computing device, alocation affected by an emergency event; generating, via the at leastone computing device, an emergency message corresponding to theemergency event; and transmitting, via a power modulator, the emergencymessage via a guided surface wave launched by a guided surface waveguideprobe.
 14. The method of claim 13, further comprising embedding, via thepower modulator, a data signal into a power signal of the guided surfacewave.
 15. The method of claim 13, further comprising exciting, via thepower modulator, a charge terminal of the guided surface waveguide probeto generate an electromagnetic field that synthesizes a wave frontincident at a complex Brewster angle of incidence of a terrestrialmedium to generate the guided surface wave that travels across theterrestrial medium.
 16. The method of claim 15, further comprisingproviding, via a feed network electrically coupled between the chargeterminal and the power modulator, a phase delay that matches a wave tiltangle associated with the complex Brewster angle of incidence.
 17. Themethod of claim 13, further comprising: receiving, via the at least onecomputing device, a message indicating the emergency event occurred froma third party service; and authenticating, via the at least onecomputing device, the third party service based at least in part on adigital signature.
 18. The method of claim 13, further comprisingdetermining, via the at least one computing device, an affected area forthe emergency event, wherein the emergency message comprises coordinatescorresponding to the affected area.
 19. The method of claim 13, whereinthe at least one computing device is further configured to at least:determine an affected area for the emergency event; and determine aplurality of codes corresponding to the affected area; generate arespective score for each of the plurality of codes; and minimize aquantity of the plurality of codes based at least on part on therespective score for each of the plurality of codes, wherein theemergency message comprises a minimized subset of the plurality ofcodes.
 20. The method of claim 13, wherein the at least one computingdevice is further configured to at least: determine an affected area forthe emergency event; and determine a plurality of codes corresponding tothe affected area; and generate at least one category code based atleast in part on the plurality of codes, wherein the emergency messagecomprises the at least one category code.